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91. (CONCEPTO DE CAMPO. CAMPO ESCALAR Y VECTORIAL. 


En general entenderemos por Campo un conjunto de valores unívocamente ligados a cada posición en un 
espacio geométrico ordinario. Se trata de una relación unívoca porque a cada punto le corresponderá en cada 
instante un único valor, aunque dicho valor podrá variar o permanecer estacionario a lo largo del tiempo, sin 
embargo no es biunívoca ya que un mismo valor podrá corresponder en cada instante a más de una posición. 


Cuando la magnitud física a la que pueda asociarse un campo sea una magnitud escalar, por ejemplo la 
temperatura, el campo correspondiente se denominará CAMPO ESCALAR y vendrá expresado analíticamente 
mediante una función del vector de posición F y del tiempo f en su caso (función escalar de variable vectorial). Así, 
si una habitación está calentada desigualmente, a cada punto de su interior corresponderá, en general, una 
temperatura diferente según su localización que podrá variar además en el tiempo. La descripción matemática de 
dicha situación corresponderá a establecer la función que haga corresponder a cada punto, en cada instante, el valor 
correspondiente de la temperatura. 


Tendremos, pues, definido un CAMPO ESCALAR en un espacio geométrico 
ordinario cuando a cada punto, en cada instante, le corresponda un escalar 
U(x yz) determinado analíticamente mediante una función del parámetro tiempo y de 


su localización geométrica en el espacio de referencia: 
== y U = U(F, £) (9.1) 


que se podrá expresar también utilizando coordenadas cartesianas como una 


Figura 9.1 E ] 
función de cuatro variables U = U(x, y,z, £). 


Análogamente, cuando la magnitud a la que pueda asociarse un campo sea una magnitud vectorial, por 
ejemplo la velocidad de las partículas de un sistema material rígido o de un fluido que diferirán, por lo general, de 
un punto a otro, el campo correspondiente se denominará CAMPO VECTORIAL y cada una de las componentes de los 
vectores del campo vendrán expresadas analíticamente mediante una función del vector de posición F y del tiempo 
í en su caso, 


Su expresión analítica se obtendrá, en este caso, mediante una función vectorial 
(de variable vectorial; es decir, de las coordenadas del punto y del tiempo): 


a - 


EA D=V(F, 1) (9.2) 
Uy2) 


o y donde cada una de las componentes del vector será una función escalar 
8 dependiente de las coordenadas del punto y del tiempo t, que en notación 
Figura 9.2 cartesiana adoptará la forma: Y = V(x, y, z, £)¡ + Vx, y,z, Oj + Ve, y,z, DK 


Tendremos, entonces, definido un CAMPO VECTORIAL en un espacio geométrico ordinario cuando a 
cada punto, en cada instante, le corresponda un vector determinado analíticamente mediante una función 
vectorial del parámetro tiempo y de su localización geométrica en el espacio de referencia. 
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Utilizando la notación compacta podríamos escribir: U = U(x,,x,,x3,1 ) y Y = V(x,,x,,x3,£) o, aún 
más sencillamente: U = U(x;, h, V = U, (x;, 1) ,, pero la más significativa puede ser la notación vectorial 
empleada en (9.1) y (9.2) donde se indica que la función depende del vector posición F e incluso extender estos 
conceptos al caso en que, en lugar del vector F, se considere otro vector cualquiera definido univocamente para 
cada punto del espacio. 


Los campos se clasifican en VARIABLES, en general cuando dependen de £, y en ESTACIONARIOS O 
permanentes cuando las funciones asociadas no dependen del tiempo sino únicamente de las coordenadas del punto, 
pudiéndose expresar en tal caso, limitándonos exclusivamente a escalares y vectoriales, como U = U(F, t) y 
Y = Y (7). En la práctica nos interesará lo fácil; es decir, conseguir que el parámetro tiempo no intervenga. Por 
ejemplo, en el acondicionamiento de un local para que la temperatura no varíe según el momento del día como 
ocurre en el exterior. 


Asimismo diremos que un campo será UNIFORME cuando todos los puntos del espacio geométrico 
considerado tengan asociada las mismas magnitudes, aunque puedan ser a su vez variables porque cambien en el 
tiempo para todos los puntos igual. Por ejemplo, en el cálculo de conductos de instalaciones nos interesará que el 
campo de velocidades sea uniforme (régimen laminar) para evitar turbulencias, aunque en este caso veremos que 
intervendrán otros factores en la realidad por los cuales aunque los vectores sean paralelos, las velocidades de las 
partículas próximas a las paredes de los conductos serán menores a las del centro. En el caso del campo de 
intensidades gravitatorias admitiremos que es uniforme sólo para ámbitos muy pequeños en los que no se aprecie 
la curvatura de la tierra (para que los vectores sean "casi” paralelos) y para altitudes muy próximas para que no 
se aprecien las diferencias en el valor de g=9,8 Nw/Kg? (mal conocida como "aceleración de la gravedad" g=9,8 m/sg?). 


92 REPRESENTACIONES GRÁFICAS: PROBLEMA REAL. 


En general, supondremos que todas las funciones que intervienen en los campos son continuas, derivables 
y uniformes, o lo que es lo mismo, le corresponde un sólo valor de la función a cada punto del espacio (obviamente 
es el punto geométrico, no el material el que tiene asociado el campo), aunque la relación no será biunivoca 
pudiendo tener varios puntos un mismo valor. 


Conocidas las funciones tendríamos perfectamente definido el campo, pero el problema real es que por 
lo general se van a desconocer dichas funciones. Se sabe que existen, pero no se conocen las expresiones que las 
definen y se llegará a ellas por aproximación a la realidad mediante interpolación de algunos valores obtenidos 
empíricamente. Obviamente su aproximación dependerá de algunos factores como el grado de exactitud en los 
cálculos, el n* de puntos (nube) medidos...etc. Hoy día, con ayuda de los medios electrónicos e informáticos de que 
disponemos lo mejor es pues guiarnos de unas INDICATRICES que nos permitan determinar, por ejemplo, cuántos 
tienen el mismo valor a un punto dado o cómo varían los valores en las proximidades de un punto. 


Para expresar de un modo gráfico un campo escalar se suele representar el lugar 
geométrico de los puntos del espacio en los que el campo adopta un mismo valor; 
es decir las llamadas SUPERFICIES EQUIESCALARES (/1g.9.3) que son las que 
resultan al dar a la función U valores numéricos definidos: 


U (x,, £) = O, 
U (x;, 1) = U, (9.5) 
Figura 9.3 U (X;, D) = U, 
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Por ejemplo, si tenemos un campo 7 = k(x?+ y? + 2? + ct?) las superficies equiescalares se obtendrán 
dando valores a 7' en cada instante / de forma que la expresión x?+ y?+ 7? = z - ct? responderá a esferas de radio 
1? = E ct?. 
k a 

Frecuentemente, las diferencias U,-U,, U,-U,, ... se procura tengan valores constantes. En el caso particular 
de que estas diferencias valgan la unidad, las superficies equiescalares corresponderían a valores de la magnitud 
escalar Up, Ugs1, Uoso» --- Up, que se mantendrían constantes sobre cada superficie. Si, por ejemplo, se tratase de 
un campo de temperatura, el dibujo de estas superficies, que se llamarían "isotermas", nos permitiría ver, de una 
ojeada, la posición de las regiones cálidas o frías. Pero su aplicación se extiende también a funciones derivadas de 
otras vectoriales o tensoriales de las que solo interviene o interesa la componente normal a un plano fijo, como es 
el caso de la presión (las isobaras nos dan los puntos de igual presión atmosférica) o de la iluminancia en el plano 
de trabajo (las curvas isolux nos dan los puntos con igual nivel de iluminación) ... etc. ya que la representación del 


campo original no siempre tiene interés o aplicación. 


Sin embargo para ciertas magnitudes es útil la representación del campo vectorial y como el módulo del 
vector del campo solamente, por lo general, no dice nada puesto que se trata de magnitudes "dirigidas", resulta 
preciso buscar otro tipo de representación. El visualizar los puntos con igual módulo de vector si éstos tienen 
dirección arbitraria no parece tener mucha aplicación pero sí podría tenerlo el saber los puntos de un entorno cuyos 
vectores tuviesen la misma dirección, aún teniendo distinto módulo. Nos aprovecharemos, por tanto, de la 
propiedad ya utilizada anteriormente de la tangencia de vectores a líneas (véase las trayectorias) para definir las 
indicatrices de los campos vectoriales como la familia de líneas envolventes de los vectores del campo o, dicho de 
otro modo, la familia de líneas a las que son tangentes los vectores del campo en cada instante, que denominaremos 
LÍNEAS DE CAMPO o LÍNEAS DE VECTOR. Para su representación (/2.9.4) trazaremos algunas de éstas líneas según 
los siguientes convenios de dibujo: 


Partiendo de un punto P, nos trasladaremos a otro muy próximo, P,, en la dirección 

S Ñ y sentido del campo en P,, Do, trazando el segmento PP, . El valor del campo en 

Y Ya  P,seráalgo diferente que en P, y describiremos un nuevo segmento P,P,, pequeño, 

a Y ?, E Y en la dirección del campo en P,, A Representado el campo en P,, Pos 
” continuaríamos el proceso, obteniendo una línea quebrada P¿P,P,...P,. El límite de 
Ús esta línea, cuando los puntos P, P, ... P, tienden a aproximarse indefinidamente 
P, entre sí, es lo que se denomina línea de vector y tiene la propiedad de que la 
dirección de la tangente en cada uno de sus puntos coincide con la dirección del 


campo en el mismo. 


Figura 9.4 


En primer lugar sabemos que son líneas porque planteamos el problema inverso: dadas las tangentes 
encontrar una curva envolvente y su ecuación matemática se obtendrá expresando el paralelismo entre el segmento 
infinitesimal de la misma drdx, ii, y el vector campo v, en el punto considerado; esto es, A dr = Y: 


Ahdx = U, (x,y,z,t) 
Ay = U,(x,y,z,1) => eliminando A => 


mI 

E 

> : y á dz = U,(x,y,z,t) (9.6) 
dl ts 


U,(x,y,2,)  U,(x,y,2,1)  U,(x,y,2,t) 
Figura 9.5 


equivalente a un sistema de dos ecuaciones diferenciales que integradas nos dan dos familias de superficies (con 
dos ctes. de integración) nd 5 és ba , cuyas intersecciones para cada valor de las ctes. y de t serán líneas que no 
podrán cortarse entre sí porque al cambiar los parámetros cambiarán las superficies pero serán superficies de la 


misma familia. 


Curso 1997/98 - Prof" I. Más - E.T.S.A.M. 


También puede justificarse porque si se cortasen en un punto existirían dos 
valores de la función para ese punto, lo cual no es posible. Otra cosa es que se 
O trate de un punto donde concurran todas las líneas de vector que se denominarán 
FUENTES O SUMIDEROS y evidentemente éstos son puntos particulares del 


campo (/12.9.6). 
o Claro que la integración analítica no tiene por qué ser nada fácil pero con los 
medios actuales se puede recurrir a la resolución gráfica por puntos que puede 
Figura 9.6 calcular el ordenador. 


Por ejemplo, dado el campo vectorial Y=(cy-bz)+(az-cx)j+(bx-ay)k haremos 
E a 0 . _€ , 293 » , 

cy=b2 az-ex  bx-ay Ue mediante el recurso de multiplicar cada numerador y denominador por el mismo valor 
nos proporcionará, al sumar todos los numeradores y todos los denominadores, unas expresiones diferenciales de 


inmediata integración: 


adx  _ _bdy  _ cdz  _ adx+bdy+cdz 
acy -abz abz -bex bex-acy 0 ax +by +CZ=C; 
de donde > 2 2 
xdx_ _  ydy  _ idz  _ xdx+ydy+idz Xx + +Z O =Ca 
CXy -bxz  ayYz1-Cxy bxz -ayz 0 


que representan la intersección de la familia de esferas de centro el origen y radio 1”=c, con la familia de planos 
paralelos dados por c, 


Por último, para indicar el sentido del campo utilizaremos cabezas de flecha 
sobre las líneas de vector y, para dar una idea de su intensidad, convendremos 
en dibujar (fg.9.7) tantas líneas a través de la unidad de área ortogonal al campo 
Es en cada punto como unidades tenga el módulo del campo en él. Naturalmente, 
P esto sólo puede hacerse de modo aproximado ya que el campo varía de un punto 
a otro sobre la unidad de superficie. De esta forma conoceremos los lugares 
donde el campo es más intenso (p.ej. P,) según las líneas de vector estén más o 

Figura 9.7 menos próximas entre sí. 


os GRADIENTE DE UN CAMPO ESCALAR. 


Dado que por cada punto pasa una superficie equiescalar, observamos que 
moviéndonos por una de esas superficies el valor del campo no varía pero si 
cambiamos de superficie obviamente cambia dicho valor. Se trata de averiguar ahora 
cómo varía el valor del campo al cambiar la posición geométrica y con qué 
intensidad lo hace en función de la dirección del desplazamiento que vaya a efectuar 
ya que, en una primera aproximación, intuimos que dicha variación dependerá de 
la posición en el espacio geométrico, es decir, nos encontramos ante un campo 
derivado del campo escalar, pero, fundamentalmente, dependerá de la orientación 
y, en consecuencia, no podrá ser una magnitud escalar. 


Figura 9.8 


Efectuando un desplazamiento en el entorno del punto (por supuesto pequeño para no perder de vista lo 
que pasa en el proceso) que llamaremos dF, admitiremos que la variación de la función escalar U es infinitesimal 
AU = dU + 5 d?U +.. y despreciando los infinitésimos de orden superior vendrá dada por su diferencial total: 
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_ QU 
— de = — dx, 
qe de Ox, ' vd 


AU= qu = Zdss Ba E 


que interpretaremos como "la suma de lo que varía en cada una de las tres direcciones proporcionalmente al 
desplazamiento dado en ellas" (variación por unidad de longitud por la componente del desplazamiento en esa 
dirección), admitiendo la superposición lineal de efectos con la arbitrariedad que ello implica. Pero no cabe duda 
que esto es solo una aproximación ya que habría que saber el valor del campo escalar en un punto para obtener el 
valor de la función en los puntos cercanos y así sucesivamente para conocerlo en todos los puntos. Para tantear su 
significado consideraremos el resultado, 


Como se trata de relacionar la variación de U con el desplazamiento diferencial dado dr y las componentes 
de este vector aparecen reflejadas en la expresión (9.7), podemos interpretarla como el producto escalar de 
di = dxi +dyj +dzk por otro vector que nos inventamos, definido por: 


9U; + ¿00 0U7 
o o d U 
rá ay) * Oz gra 12) 


que denominamos gradiente del campo por indicar una variación derivada de un desplazamiento, ya que resulta: 


dU = gradU: di (9.9) 


$e agro gradiente podrá expresarse también introduciendo un vector simbólico (nabla) 
y = Z ¿+2 7+2É enla forma VU= e il, = gFadU = Y y dado que U es función de las coordenadas del punto 
(función escalar de variable vectorial), sus derivadas parciales y, por tanto el vector gradiente, también lo serán, 
obteniéndose, en general, una función vectorial de variable vectorial. Es decir, si bien obtendremos, particularizado 
a cada posición, el vector gradiente del campo en el punto considerado, acabamos de determinar un campo vectorial 


especial por derivar de otro escalar, 


Obviamente dU será nulo si nos movemos sobre la propia superficie equiescalar 
que pasa por el punto para el que particularizamos el gradiente, ya que U (fig.9.9) 
tendrá el mismo valor. En tal caso dr, será un vector contenido en el plano 
tangente a la superficie por dicho punto y, dado que, en general, el gradiente no 
será nulo, la dirección de gFadU tendrá que ser ortogonal al plano para que sea 
nulo el producto escalar de ambos vectores dU = gradU* dF, = 0 coincidiendo, 
por tanto, con la normal a la superficie equiescalar en cada punto. En 
consecuencia, las líneas de vector del campo de gradientes serán también 
ortogonales a las superficies equiescalares constituyendo lo que en matemáticas 
se denomina “trayectorias ortogonales a la familia de superficies U”. 


Figura 9.9 


Para conocer el sentido de gFad U nos moveremos en la propia dirección del vector gradiente y como es 
arbitrario diremos que su valor será positivo cuando tenga el sentido de los valores crecientes de U ya que, 
denominando dn al módulo de 47, podremos escribir: 


dU, = | gradU| |d,|cos |? 


. z =+ | gFad U| dn 0.10) 


donde el producto escalar se ha reducido al producto de los módulos por tratarse de vectores con la misma 
dirección. 
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Por último, para averiguar el módulo de g*ad U buscaremos igualmente el avance en la propia dirección 


del vector para eludir el parámetro en la expresión |2rad U | = pa por lo que tendremos: 
P|cosp 
ms dU 
leradU| = — (9.11) 
dn 


que podríamos interpretar como la variación del campo U “por unidad de longitud” según la dirección de la 
normal a la superficie equiescalar en cada posición. Por tanto, si introducimos el vector unitario en dicha dirección 
y sentido el de los valores crecientes de U podremos escribir la expresión vectorial del gradiente en la forma: 


gradu= Ly di 
dn 


n 


cuyo significado será el de un vector que, en cada punto donde la función U y sus derivadas primeras estén 
definidas, indicará la dirección y sentido de máxima variación del campo escalar, así como su valor. 


Así, por ejemplo, dado el campo escalar U =x?+y?+2?, estacionario al no depender de t, cuyas 
superficies equiescalares U = 1? son las esferas con centro el origen de referencia, el campo de gradientes será: 
gradU =2xi+2yj +2zk = 27 es decir, el vector gradiente en cada punto tendrá la misma dirección del vector de 
posición correspondiente (de donde se deduce que las líneas de vector del campo serán el haz de rectas que se 
cortan en el propio origen de coordenadas) y su módulo será 2r como se comprueba al ser - =2r =|grad U| . 


No obstante la expresión (9.9) permite también determinar cuánto varía U en 
cualquier otra dirección ya que realizando un desplazamiento direccional 
infinitesimal dr que forme un ángulo o. con la dirección del gradiente será: 
dU = gradU-: cosa: dr, de donde obtenemos la “variación por unidad de 
longitud” en esa dirección o, más adecuadamente, la derivada direccional: 


dU 

— = gradU:coso. = grad U = gradU: € 

Figura 9.10 dr g gra0, £ % Py 
expresable también simbólicamente como grad , U que se denomina "gradiente de U según la dirección que forma 
un ángulo a. con la de máxima variación", que equivale a proyectar gFad U sobre la dirección de dr (fig.9.10) o, 


lo que es igual, a multiplicar escalarmente el vector gradiente por el unitario correspondiente. 


Si se considera como magnitud vectorial, la derivada direccional se deberá multiplicar por un vector 
unitario según la dirección considerada y el lugar geométrico de estos vectores será una esfera tangente a la 
superficie equiescalar cuyo diámetro es grad U, observándose que la definición (9.1/) inicial de gradiente 
corresponde a sumar vectorialmente el valor del gradiente de U según cada una de las tres direcciones que 
corresponden a los ejes. Si uno de los ejes, p. ej. el ox, es ortogonal en un punto a alguna de las superficies 
equiescalares, el gradiente en dicho punto será — ¡, expresión equivalente a la (9.12). Esta última fórmula está 
matemáticamente bien definida si se define U en función de las distancias al punto considerado según la normal 
y dos direcciones ortogonales entre sí situadas en el plano tangente. 


Nótese, por último, que el resultado sería totalmente distinto si tratásemos de averiguar la variación de 
un campo vectorial AY = dY + > d?V +... Despreciando igualmente infinitésimos de orden superior tendríamos: 


CIAO PA 


A ay de 

AV < dv dev 7+v,k] A A > 3 grad V> dr 
9, av, or, | dz 
a y e 


es decir, nuestro "vector" gradiente tendría 9 componentes siendo necesario recurrir a la notación tensorial. 
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94  ELUJO DE UN CAMPO VECTORIAL. DIVERGENCIA. 


De modo análogo a cómo el campo de gradientes de un campo escalar relaciona, mediante una expresión 
analítica, la distribución espacial del campo según superficies equiescalares con la variación de la función, a 
continuación analizaremos la correspondencia entre un campo vectorial y otro campo (en este caso escalar) que 
permitirá relacionar analíticamente la distribución espacial de aquél según las líneas de vector con la variación del 
campo vectorial en el entorno de cada posición. 


Para ello, dado un elemento de superficie AA representado por el vector normal 
cuyo módulo coincide con el área y su sentido se relaciona con el de recorrido del 
contorno de la superficie, según la regla del sacacorchos (7ig.9.1/), necesitamos 
definir el "/lujo elemental" de un campo vectorial Y a través del elemento de área 
AA, como el resultado de realizar el producto escalar: 


Ap = V- Mm (9.15) 
Figura 9.11 


donde Y es el valor de la función vectorial en un punto de AA, que puede ser cualquiera pues el elemento de 
superficie se supone muy pequeño (en el límite, infinitesimal, dA). 


Si se trata de una superficie finita A, la dividiremos en pequeños elementos de 
área AA, (fig.9.12) y, considerados como vectores, obtendremos el flujo a través de 
A como límite de la suma de los flujos elementales a través de todos los elementos 
de superficie; es decir: 


n 


9 = lim 2 PV, M4, = [[ V> di (9.16) 
A 


Figura 9.12 pl 

Como se indica, el límite no es otra cosa que la integral sobre la superficie (integral doble) del producto 
escalar Y - 4 y generalmente es finito pues, aunque el tamaño de cada elemento tiende a cero, su número aumenta 
indefinidamente. 


Si se ha dibujado el campo vectorial de acuerdo con las prescripciones indicadas 
en $9.2, el resultado de calcular el flujo a través de una superficie coincide con 
el número de líneas de fuerza que la atraviesan. Efectivamente (/18.9.13), el 
producto escalar puede escribirse: 


V-M=V-:M:-c0s0=V-: A, (9.17) 


Figura 9.13 =, 
donde AA, es la proyección de AA sobre el plano normal al campo Y. 

Como el módulo de Y mide el número de líneas de campo que atraviesan la unidad de área transversal, 
el valor VAA, será la contribución del elemento AA al número de líneas que atraviesan un área finita y la suma 
(9.15) nos da el número total de líneas que atraviesan la superficie A, En este sentido cabe definir lo que se 
denomina "tubo de corriente" como la superficie formada por el conjunto de todas las líneas de vector que pasan 
por los puntos de una línea cerrada. Si esta línea es sencilla, dicha superficie tiene, en efecto, la forma de un tubo. 
Evidentemente, el flujo a través de cualquier elemento de superficie de uno de estos tubos es nulo, pues A4 ¿ Será 
ortogonal al vector del campo Y. 
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= o - En el caso de una superficie cerrada, donde se toma como positivo el sentido de 
la normal exterior, el ángulo o, puede ser superior a 90, en cuyo caso el flujo se 
contaría negativamente. Se entenderán, por tanto, como positivas las líneas 
"salientes" (flujo saliente) que en la fig. 9.14 pasan a través del área FHG y como 
negativas las "entrantes" (flujo entrante) que atraviesan el área FGL, de manera 
que al calcular el flujo total a través de la superficie obtendremos la diferencia 
entre el número de líneas que salen y que entran, es decir, las que nacen en el 
interior de la superficie si el flujo total es positivo (o mueren si es negativo). 


Figura 9.14 


En la figura 9.14, el flujo será positivo e igual aproximadamente a 2, pero en otra posición el flujo del 
campo a través de la misma superficie podría adoptar otro valor en función de sus características, por lo que parece 
útil definir una magnitud que establezca una relación entre el número de líneas de vector que nacen o mueren en 
el entorno de cada posición, independientemente de la superficie a considerar. 


Para obtenerla analíticamente tomaremos una superficie pequeña encerrando un 
volumen At y dividiremos su área en un número » relativamente alto de 
elementos de área AA; (fig. 9.15), de forma que la expresión » dt tendrá un 
valor, en general, finito aún cuando supongamos que la subelficie se contraiga 
en el entorno de un punto genérico P de forma que At -> 0, 


Figura 9.15 
Como el numerador mide la diferencia entre el n* de líneas de campo 


que salen y entran de la superficie y dividimos por el volumen encerrado en ella, dicho cálculo correspondería al 
número de líneas de campo que nacerían (o morirían) por unidad de volumen si el espacio encerrado fuera lo 
suficientemente pequeño para admitir que apenas variaría el valor del campo en él. Por tanto, cuando nos limitamos 
al entorno del punto haciendo At - 0, el valor obtenido será una característica del campo en ese punto que no puede 
depender de la forma de la pequeña superficie considerada, ni de cómo se realice la operación de límite en torno 
aP, obteniéndose un resultado unívoco que se denomina DIVERGENCIA del campo vectorial en el entorno de P y 
se expresa como: 


e ” Vo-A, 
div Y = lim Y) —— (9.18) 
Ar-0 1 At 


Es decir, si el resultado del cálculo indicado fuese, por ejemplo, a, esto significaría que si en la unidad de volumen 
el campo tuviese las características que tiene en un entorno de P, nos sería preciso hacer nacer "a" líneas por unidad 
de volumen para dibujarlo adecuadamente. Si no nacen ni mueren líneas de campo en una región del espacio, en 
ella se verificará que el flujo total a través de cualquier superficie cerrada será nulo // Y- dA = 0 pero si el campo 
vectorial del que procede es de tal naturaleza que no se precisa, para representarlo, hacer nacer (o morir) líneas de 
campo, por prolongarse indefinidamente o por cerrarse sobre sí mismas, se verificará que en todos los puntos 
div Y = 0 y ésta será la característica matemática de un campo SOLENOIDAL. 


La magnitud que se obtiene es, evidentemente, un escalar asociado al punto P que puede variar de un punto 
a otro del espacio en función de cómo varíe el campo vectorial y así como la operación gradiente nos permite 
asociar un campo vectorial a uno escalar, esta operación, generalizada a todos los puntos del espacio a que se 
extiende el campo vectorial Y, nos permite asociar un campo escalar a otro vectorial ya que, en general: 


divV(F,) = 0(F,1) (9.19) 


aún cuando pueda obtenerse un campo uniforme (independiente de la posición) y permanezca estacionario a su vez 
en el transcurso del tiempo como veíamos para todo campo escalar. 
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Para encontrar la expresión de la divergencia en coordenadas cartesianas y dado que el resultado de aplicar 
la expresión (9.18) no depende de la forma de la superficie, supongamos que se trata de un sencillo paralelepípedo 
cuyas caras son paralelas a los planos coordenados y que Y es el valor del campo en el punto central. 


A través de la pareja de caras paralelas al plano 0x,x3 solo producirá flujo la 
componente V, ortogonal a ellas (.9.16), cuyo valor, para la cara anterior y 
posterior respectivamente, será, despreciando términos en Ax: 

WV, Ax; 7% VW, Ax; 


V, + —— a 
E: ra 2 


resultando para el flujo de la componente V,: 


Figura 9.16 


Al considerar las tres parejas de caras paralelas tendremos para expresión del flujo total: 


Ap = EV, M4, = (E ' E ' =] Ax, Ax, Ax, y como el volumen del paralelepípedo es Atr=Ax,Ax,Ax,, quedará: 
2% a 
=$ ov, 9V, oV, 9v, 
dVe As Ll4 242 2 E O E (1,3, + V,0,+ 1,12) = 9-9 
al dx, dx 0%  0x a” dE da > a dd E (9.20) 


que es la expresión buscada mediante el operador nabla definido con anterioridad. 


9s CIRCULACIÓN DE UN CAMPO VECTORIAL. —ROTACIONAL. 


Para una línea AB cualquiera, situada en un campo vectorial, el vector del campo Y adoptará, por lo 
general, un valor diferente en cada uno de sus puntos, por lo que descomponiéndola en elementos de arco muy 
pequeños de forma que pueda admitirse despreciable la variación de Y en cada uno de ellos (/12.9.17) definiremos 
circulación elemental del vector como el producto escalar AC = 7 -AF. 


La suma de las circulaciones elementales para todos los elementos de arco definidos entre A y B se 
denomina "circulación del vector Y a lo largo de la línea AB" y se expresará en el límite, cuando AF- 0, 
mediante la integral: 


B 
C=limYEV +47, = [Vd (9.21) 
Ár¡-0 A 


o bien, si y? está expresado en cartesianas Y = Y, (x)3, + V,(c)4,+ Vx 3, = V, 1, 


B B 
C = f V,dx, + V,dx,+V,dx,= f V,dx, (9.22) 

Figura 9.17 A 4: 
En el caso de que el integrando pueda expresarse como diferencial de una función escalar U, esta integral 
no dependerá del camino recorrido, sino únicamente de los puntos A y B inicial y final del mismo, ya que entonces: 


B B 


> 


B 
C= fV:dr= fdU=Ug=- U, (9.23) 
a A A 


pudiéndose considerar el campo vectorial como campo de gradientes de un campo escalar ya que, de acuerdo con 
la (9.9), V+ dF = gFad U* di = dU. 
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Diremos entonces que Y*= gFad U es "conservativo" pues si la línea es cerrada 
(fig.9.18) tendremos: 


C=4V"-dP=U,-U,=0 


A No ocurrirá igual si la función no tiene primitiva (ni siquiera es integrable) 
dr debiendo sustituir las variables en función de los parámetros de la curva descrita 
Figura 9.18 que serán distintos cuando el punto pasa por A y por B=A. 


Mediante un proceso análogo al seguido para la definición de la divergencia a partir del flujo de un campo 
vectorial intentaremos definir ahora una nueva magnitud que denominaremos ROTACIONAL del campo vectorial 
para expresar la circulación del vector del campo en el entorno del punto, para lo cual necesitamos delimitarlo 
mediante unas superficies elementales arbitrarias que en el límite, cuando A4A,- 0, encerrarán el propio punto P. 


Dado que la circulación del vector del campo a lo largo de la línea que encierran una de esas superficies 
elementales es un producto escalar, el valor que obtendremos “por unidad de superficie” = no resultará 
representativo de su orientación a pesar de que resulta obvio que así será, por lo que intuimos ya que dicha 


magnitud habrá de ser vectorial. 


Consideraremos, por tanto, una línea plana y cerrada que limite un área A4 
situada en un campo vectorial Y y tomando un sentido de recorrido la 
dividiremos en pequeños segmentos de cuerda AF,, calculando el producto 
V.: AF, que, en el límite, será la circulación (distinta de cero si el campo no es 
conservativo) del vector Y a lo largo de dicha línea. Definiremos la componente 
sobre la normal al plano (con el sentido que corresponda aplicando la regla del 
sacacorchos a Ar,) del nuevo vector que llamaremos Rotacional de V a la 
cantidad: 


Figura 9.19 rot, Y = lim = lim L = rotV-5¡ (9.24) 


que puede entenderse asimismo como la proyección de ese vector denominado rotacional sobre la normal o, en otras 
palabras, como el producto escalar por el unitario según la normal. 


Si P es un punto en cuyo entorno se mantiene la línea para AA - 0, podremos trazar otros dos planos 
-de modo que P sea el vértice de un triedro trirrectangular- y trazando sendas líneas cerradas sobre ellos 
calcularemos, mediante (9.24), las componentes, según las direcciones ortogonales a ambos, de rót Y, obteniendo 
así las tres componentes, según las aristas del triedro, que pudieran definir un vector. Esto no es evidente ya que 
se requiere que dicho vector resulte el mismo aunque se le defina a partir de otros planos cualesquiera de los 
infinitos que pueden pasar por P. 


Para comprobar que así ocurre tomaremos (fig.9.20) los triángulos PRS, PSQ, 
PQR como las líneas que tienden a reducirse al entorno de P, con lo que se forma 
un cuarto triángulo QRS que, al disminuir los otros tres, también quedará en el 
entorno de P. Sean ñ,, ñ,, ñ, y ¡ vectores unitarios normales a los planos de los 
triángulos considerados y C,, C,, C, y C las circulaciones a lo largo de los 
triángulos PRS, PSQ, PQR y QRS, respectivamente, con lo que el rotacional 
podrá escribirse, a partir de los tres primeros triángulos: 


e Ms ÓN 
Figura 9.20 rót V = lim y a ds (9.25) 
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Comprobaremos que éste vector tiene como componente, según la normal X al triángulo QRS, el valor de dicha 
componente tal como se obtiene considerando la definición (9.24) aplicada a dicho triángulo. Efectivamente, 
multiplicando escalarmente (9.24) por 7, tendremos: 


rot. D = liml L cosa, + cosa, + cos a 
no 7 TA 157, 27, 3 (9.26) 


donde a,, a,, a, son los ángulos que forma ; con las normales a los restantes triángulos, que coinciden con los 
ángulos de los planos que los contienen, de forma que, si A4 es el área del triángulo QRS, se verifica: 
M4 = 2 ._% ._% pudiéndose escribir la (9,25) en la forma: rot,7 = lim £0%:9 ño 3. Ahora bien, en la figura 9.20 
se observa que, si recorremos los triángulos del triedro de forma que SQR sesilte recorrido positivamente, los 
términos de circulación correspondientes a las aristas PQ, PR, PS se recorren dos veces en sentido opuesto, luego 


C¡+C,+C,= C, obteniéndose: ror,? = timÉ. que coincide con la definición (9.24) aplicada al QRS. 


se En consecuencia, rot V se comporta, efectivamente, como un vector cuya 
Ma v : ' a Ñ ' 
expresión en cartesianas podremos deducirla sin más que considerar la 
pen circulación de Y según una línea de forma sencilla situada en (/ig.9.2/) cada uno 
de los tres planos coordenados (p. ej, el rectángulo de lados x,, x, para el plano 
ox,x, nos dará la componente según 0x,). Así, si Y es el vector del campo en el 


0) j punto central P, obtendremos para la circulación C,, en una aproximación lineal 
para las Ax; 


4 


Figura 9.21 


Y, Á WV, A wm, wV, 4 o, vw, 
MS 


por lo que, teniendo en cuenta que el área del rectángulo es AA=Ax,Axy y aplicando la (9.24), se obtiene para 
componente del rotacional según 0x;: 


Si consideramos ahora otros dos rectángulos centrados en P y en planos paralelos a los 0x,X; y 0X,X), 
Y  %W 
obtendremos las otras dos componentes sin más que permutar los subíndices circularmente: — rol, fa - 


A h %, % 
rot, V = E Ss en resultando para expresión final del rotacional: 
1 
a o 0v. 9Y,  0v,. 9,  oY 
rotV = | 2 - lg, +| 2 - —l8,+| — - 18, (9.27) 
OX, Ox, Ox,  0x, Ox, Ox, 
que puede escribirse como un producto vectorial simbólico: 
ZE KE E 
e 0 08 09 a 
0 <= <A ER ls RP (9.28) 
Ox, 0X, Ox, 
Y, v, V, 


sil» 
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9.Ó 


[m] 


(0) 


E 


CONCEPTO MATEMÁTICO DE GRADIENTE, DIVERGENCIA Y ROTACIONAL.. 


- Si, como hemos visto, TODO CAMPO ESCALAR U TIENE ASOCIADO UN CAMPO VECTORIAL "derivado" 


de él (y no le llamamos así por casualidad): 
U(xyz) > gradiente > V=gradU=VU 


nos surge la pregunta si se cumple también a la inversa; es decir, si todo campo vectorial es asociado de 
un campo escalar o, más exactamente, si todo campo vectorial deriva de un campo escalar. Una forma 
inmediata de comprobarlo es observar qué sucede con las componentes del campo vectorial ya que si 
efectuamos las segundas derivadas de las componentes de vector gradiente se igualan las derivadas 
cruzadas al no influir el orden de derivación de una función escalar, y al ser: 


grad U = + EA =k = V(xyz)í+ V(xyz)j + V(xyz)k =P 


tendrá que cumplirse también la igualdad de las derivadas cruzadas entre las componentes del campo: 


$0 _ 40 a al -2(2) a BA 
Rx yo Y Y xXx 9 Y da 
Teorema 
7 2 y 218) ¿(2 s Us ZR 
Schwarz RX Xd % Pa óx %% % o 
$0 80 2 q a 10% E 
AD Do Y % ó% Dd DY A 


podemos afirmar que NO TODOS LOS CAMPOS VECTORIALES DERIVAN DE UN CAMPO ESCALAR. 


En consecuencia, distinguiremos a unos y otros denominando CAMPOS CONSERVATIVOS a los campos 
vectoriales que procedan de un campo escalar por aplicación del "operador" gradiente y CAMPOS 
DISIPATIVOS' a todos los demás (que, por supuesto serán muchos más). 


Ya que no podemos afirmar que todo campo vectorial derive de un campo escalar, trataremos de averiguar 
si podemos asociar un campo escalar a otro vectorial. Haciendo uso del operador nabla 
V = Z ¡+ > j+ = k, definido anteriormente para expresar el gradiente en forma abreviada -que no 
es otra cosa que un vector simbólico que opera como tal con lo que tiene a su izquierda y actúa 
exclusivamente sobre lo que tiene a su derecha efectuando las derivadas parciales correspondientes- y dado 
que el producto escalar de dos vectores es un escalar, PODEMOS ASOCIAR A TODO CAMPO VECTORIAL 
UN CAMPO ESCALAR procedente del producto escalar de dicho operador con el vector del campo, que 
denominaremos divergencia del campo vectorial: 


El por qué de dicha denominación lo comprenderemos al estudiar el trabajo de las fuerzas y la transformación de energía mecánica 
ya que si las fuerzas (vectores) no derivan de un campo escalar "disiparán" energía por rozamiento, calor ... mientras que si podemos 


encontrar un campo escalar asociado al campo de fuerzas, transformarán exclusivamente energía potencial en cinética y viceversa. 
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F 2 divergencia > (xyz) 


de forma que 


Q = V-F = divF = jLEÉ AE-I+EJ+FKk)= ci MA (9.29) 
J E 


Como hemos visto en $9.4 los campos (obviamente los FP) para los cuales divF=0 se denominan 
Solenoidales y serán muy importantes en el estudio de las deformaciones puesto que cuando lo sean no 
se producirán cambios de volumen. 


Pero como ocurría con los conservativos, no todos los campos vectoriales serán solenoidales. 
Evidentemente no todo campo escalar deriva de un campo vectorial por aplicación de la divergencia o, 
al menos, no de uno solo, dado que podemos encontrar más de un conjunto de funciones escalares que 
sumadas nos den la que tenemos de partida y como además tendríamos que integrar los tres sumandos 
elegidos, surgirían otros tres parámetros adicionales reflejados en las constantes de integración respectivas. 
Es decir, encontraríamos infinitos campos vectoriales que lo cumplirían y por tanto no merece la pena 
estudiarlos. Lo interesante es lo particular; la notabilidad es que no haya divergencia, no que sí la haya. 


Que sea conservativo el campo vectorial no implica que sea solenoidal, aunque puede darse el caso que 
un campo vectorial sea simultáneamente conservativo y solenoidal. 


Intentaremos ahora asociar un campo vectorial a otro vectorial, aunque de todas las posibilidades 
existentes solo mencionaremos una que tiene interés por las aplicaciones que va a tener en nuestro ámbito 
de estudio. Y siguiendo con el juego matemático que nos proporciona el operador nabla PODEMOS 
ASOCIAR AL CAMPO VECTORIAL DADO OTRO CAMPO VECTORIAL resultante de multiplicar 
vectorialmente dicho operador con él, basándonos en que han aparecido magnitudes en la Física que 
tienen un comportamiento similar al que describe el vector del campo resultante que denominaremos 


rotacional: 
F 3 rotacional > W 
de forma que: 
PIE 
> + or OF > 9F, OF, > 9F, or, > 
W=-VxÉ= rot Ñ= = Zz (2-2) (2-2)7 (2-2); 
x dy 0z Oy Oz Óz Ox Ox dy 
F, F, F, 


entendiendo que las "multiplicaciones" se leen "derivadas". 


Ahora bien, dado un campo vectorial W , por lo mismo que en el caso del campo de gradientes, aunque 
la demostración es más compleja, NO siempre podemos decir que deriva de la aplicación del rotacional 
a otro campo vectorial. Sin embargo hay una particularidad que merece destacar si tenemos un campo 
vectorial 4 cuyo rotacional sea nulo rot F' = 0 en todos los puntos, en cuyo caso llamaremos irrotacional, 
y es que al ser nulas las componentes del vector rotacional tendremos: 
Me 0 E, E _4 
dy ó E E 


a 
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particularidad que ya habíamos visto con anterioridad en el caso de los Campos Conservativos, por lo que 
podemos afirmar que TODO CAMPO IRROTACIONAL ES CONSERVATIVO Y VICEVERSA. Su interés será 
patente en el estudio de fluidos ya que si el campo de velocidades de las partículas es irrotacional veremos 
que éstas se van a trasladar exclusivamente, sin producirse giros (remolinos). 


| Nos queda, por último, averiguar cómo se puede pasar de un campo escalar a otro escalar, pero ya resulta 
obvio con lo visto hasta ahora, ya que si: 


U(xyz) =D gradiente > F y P > divergencia >  Q(xyz) 


con ambas operaciones simultáneas tendremos: 


U(xyz) <= divergencia del gradiente >  Q(xyz) 


cuya expresión resultante será: 


V-(VU) = WU = div grad U = 1444, - p(ayz) (9.31) 
ar? dy 97? 


conocida como Laplaciana de U, que nos conducirá a la expresión de la transmisión del calor o a la de la 
membrana en Estructuras, 


Si div (gFad U ) = 0 significa que el campo vectorial E gFad U* es SOLENOIDAL y como deriva de una 
función original U* por aplicación del gradiente, es al mismo tiempo CONSERVATIVO, por lo que a dicha 
función original se la denomina FUNCIÓN ARMÓNICA (que descubriremos en la Función de tensión o de 


Airy). 


9.7 TEOREMA DE GAUSS ed FÓRMULA DE STOKES 


Consideremos ahora una superficie cerrada A (/1g.9.22) y dividamos el volumen que encierra en 
pequeños elementos At, a cada uno de los cuales podremos aplicar la definición (9. 18) en la forma: 
lim div V - At = E V, . A, (9.32) 
At-+0 

Está claro que, si sumamos para todos ellos en el límite, el primer miembro nos dará la integral, 
extendida a todo el volumen, de div Y. En cuanto al segundo miembro, debe observarse que, al 


sumar sobre las caras comunes a dos volúmenes elementales contiguos, para cada uno de ellos será 


normal exterior lo que para el otro es interior, anulándose dos a dos todos los términos salvo los de 


Figura 9.22 las caras que componen la superficie exterior dada que encierra el volumen total. 
El resultado pues será, en el límite, la integral del flujo del vector Y a través de la superficie A, obteniéndose, en definitiva: 


ff div V- di = $$ V-di=8Q (9.33) 


fórmula de gran importancia, al relacionar una integral doble con otra triple, que es conocida como "fórmula de Gauss". Esta fórmula podría 
admitirse sin necesidad de la demostración indicada pues equivale, conceptualmente, a igualar el número de líneas de campo que nacen en el 


volumen (primer miembro) y la diferencia entre las que salen y entran (segundo). 
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Mediante un razonamiento análogo al utilizado para la obtención de la fórmula anterior, podemos 
imaginar una superficie que se apoye en una línea cualquiera y descomponerla (fig.9.23) en 


elementos de superficie, a cada uno de los cuales resultará aplicable la expresión (9.24) en la forma: 


li > S . pan > a "A cd . > z = 
lie rot,V NA de rot Y - A4 e 2 Y * Ar, (9.34) 


donde se observa que, para una superficie y su contorno que tiendan a cero, el flujo del rotacional a 


Figura 9.23 través de la superficie coincide con la circulación del vector a lo largo de su contorno. 


Al sumar para todos los elementos de superficie, en el límite, el primer miembro se transformará en el flujo del rotacional de Y a 
través de la superficie total, /rót Y - dA; en cuanto al segundo miembro, será necesario realizar la suma para los contornos de los elementos 
de superficie, pero, como todas las líneas interiores serían recorridas dos veces en sentidos opuestos, las circulaciones correspondientes se 
anularían en la suma, quedando como resultado de ésta, únicamente, la circulación a lo largo de la línea contorno de la superficie total. 


Obteniendo así la importante "fórmula de Stokes": 
ff 7-di- $7-a (9.35) 


dA Dado que la circulación del vector Y alo largo de una línea determinada tendrá también un valor 


determinado, el flujo de su rotacional será el mismo para todas las superficies que se apoyan en ella. 

Dos cualesquiera de estas superficies constituirán una superficie cerrada (/g. 9,24), por lo que, estando 

7 dirigido el sentido positivo de la normal de acuerdo con el de recorrido de la línea, el flujo de los 

rotacionales será el mismo a través de una y otra superficie. Sin embargo, considerándolas como una 

Figura 9.24 sola superficie cerrada, el signo positivo de la normal y de los vectores dA corresponde a la normal 
exterior, por lo que debemos invertir el sentido de la normal de una de las dos superficies. 

Esto significa que el flujo del rotacional será igual y de signo opuesto para una de las dos superficies. En consecuencia, el flujo a 
través de una superficie cerrada de un vector que puede expresarse en la forma W = rót Y debe ser nulo, con lo que, apoyándonos en la fórmula 
de Gauss (9.33), podemos escribir: $$ rotV- da = i 1 S div (rot Y) dí = 0 y, Como esto debe ocurrir para cualquier superficie cerrada, tendría 
que verificarse, de modo general, que div (rot Y) = 0 es decir, el rotacional de un campo vectorial es siempre un vector solenoidal. Por otra 
parte, si el rotacional de un vector es nulo en todos los puntos del campo (campo irrotacional o conservativo), se verificará, por la (9,35), 
$ V «dr =0, lo que prueba que el integrando es la diferencial exacta de una función dU = gFad U: dr, es decir Y = gFad U, por lo que 


rót ( gFad U) = 0 y el campo Y será conservativo. 


En resumen, las propiedades de los campos vectoriales /RROTACIONALES son: 


P Su rotacional es cero: yor Y = 0 
27 Su circulación a lo largo de una curva cerrada es nula $ V-dr=0 
si Son gradiente de una función escalar: Y = grad U 


Y las de los campos vectoriales SOLENO/DALES son: 


e Su divergencia es nula: diy Y = 0 
22 Suflujo a través de una superficie cerrada es cero. $ $ V-dl =0 


ge Son el rotacional de una función vectorial: Y = rot W 


cl 
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CAMPOS TENSORIALES Y CÁLCULO UENSORIAL 


a OTRO CONCEPTO DE CAMPO 


En el estudio de determinados fenómenos, como la transmisión de acciones en un medio continuo natural 
es decir, más o menos deformable, nos vamos a encontrar con algunas magnitudes direccionales cuyos valores 
vendrán dados, entre otras variables, en función de una dirección normal al plano o superficie a considerar. Sin 
ir más lejos, para no entrar en materia de otros cursos, en la mecánica del sólido rígido tendremos que relacionar 
el momento de las fuerzas actuantes con la variación del momento cinético mediante un operador indicativo de la 
distribución de la masa en función de la dirección del eje instantáneo de rotación, 


Nos encontramos, por tanto, con otro concepto de campo de un orden superior 
al que responderán estas “funciones” cuya característica común va a ser que a 
cada punto del espacio geométrico no le va a estar asignado un único valor, como 
en los campos escalar y vectorial, sino que le podrán corresponder «> valores, de 
los que se adoptarán en cada caso solo uno en función de otra variable adicional 
que va a ser una dirección definida por un vector. Es decir, las magnitudes 
asociadas en estos campos “a cada punto” dependerán, en general, además del 
tiempo t y de la posición Y del punto como en los anteriores, de otro vector Y, 
indicativo de las posibles orientaciones a considerar en el entorno de aquél. 


AS va) 


Figura 10.1 

En otras palabras, podríamos decir que el vector Y obtenido en estas “transformaciones” tendrá unas 
componentes que diferirán, no solo en función del punto considerado en el espacio SeomÉnies ordinario, sino 
también en función de las componentes del vector Y relativo a una dirección variable, wW =w (Y, F, £), lo que 
podríamos representar, para cada posición, mediante tres funciones w, =f (v,, V,, V,), w, =£(v, ¿Vs Yo) Y 
w, =£(v,, v y» v,), de forma que si dichas funciones son lineales, es decir pee melena una combinación 
lineal de las v, por unos valores escalares Ti: 


e = TY +71) + Tv, 
E TapY] + Ta, + Ta, (10.1) 
bj y Tay + Ta 2 + Tv, 
lo podríamos reducir a una notación compacta, w, a T,,», o simplemente w, = = Tv (notación de Einstein), aunque 
también resulta útil la notación matricial, il cada componente del sectors transformado W como suma 
de los productos entre los elementos de cada fila de la matriz del tensor por los del mismo subíndice en la columna 
de las componentes del vector Y, es decir: 


ME Ze Tal 
w [=[Z), Tao Ta] [va (10.2) 


*s La, ES Ts, ES 


expresable simbólicamente mediante la notación matricial reducida |w| =[]7'| |v| o bien mediante la notación 
vectorial y - 7. yy, en la que se adopta la misma simbología del producto escalar por similitud. 
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En principio, como se trata de aproximarnos a los fenómenos reales, la mayoría de las veces de forma 
empírica, se parte de una función lineal. Bastaría, por tanto, conocer los transformados de tres vectores aplicados 
en el mismo punto para determinar las componentes del tensor y si los resultados obtenidos aplicados a un cuarto 
vector se aproximan a la realidad diremos entonces que el fenómeno que representan es un fenómeno lineal (ya 
que si los T; fueran a su vez función de v,, v, y v, la cosa se complicaría aún más y el tensor ya no sería lincal). 


En consecuencia se tendrán 3? = 9 incógnitas (las componentes T;) para definir un tensor de este tipo, 
a diferencia de las 3' = 3 incógnitas necesarias para definir un vector o de las 3% = 1 necesarias para definir un 
escalar, siendo el 3 indicativo de la base trirrectangular elegida. Podríamos decir, por tanto, que los tensores que 
estudiaremos serán operadores lineales de orden 2 que transformarán un vector en otro de dirección distinta, por 
lo general, al anterior y los representaremos con dos líneas encima de la letra representativa de la magnitud 7 de 
modo análogo a la representación vectorial Y que entenderemos como operador o tensor de orden 1 y a la 
representación escalar U, entendida a su vez, como operador o tensor de orden 0. 


El conjunto de los 3?=9 valores (escalares) T que definirán el operador en un punto P, cada uno con un 
significado particular correspondiente a la magnitud física en cuestión, lo podemos denominar TENSOR del campo 
asociado a P y se entenderá que el vector origen y el transformado por él tendrán el mismo punto de aplicación. Por 
tanto el conjunto de los tensores asociados a cada uno de los puntos del espacio definirán un CAMPO TENSORIAL 
que estará constituido por 3%=9 funciones que, por analogía con las definiciones de los otros campos, expresaremos 
como función ZENSORIAL de variable vectorial, 7 =T(+,1). 


Tendremos, por tanto, definido un CAMPO TENSORIAL en un espacio geométrico ordinario cuando a 
cada punto, en cada instante, le corresponda un tensor determinado analíticamente mediante nueve funciones 
escalares del parámetro tiempo y de su localización geométrica en el espacio de referencia, mediante el cual 
a cada vector aplicado en él le corresponderá un único vector de dirección de dirección distinta por lo general. 


Y generalizando se entenderá por CAMPO TENSORIAL de orden »n el conjunto de las 3” funciones 
necesarias para definirlo en su totalidad:. 


=] >» campo tensorial de orden O ESCALAR 
S=3 2 campo tensorial de orden 1 VECTORIAL 
$=9 >) campo tensorial de orden 2 TENSORIAL 


10.2 ÉL TENSOR DE ENTORNO CAM CAMPO. 


Como ya quedó indicado al estudiar el gradiente de un campo escalar, el resultado resulta totalmente 
distinto cuando tratamos de averiguar la variación de un campo vectorial en el entorno de los puntos debido a que 
cada una de las componentes del vector de campo es una función escalar, 


Recordaremos que a todo campo escalar se le podía asociar otro vectorial denominado campo de gradientes 
ya que a cada punto del espacio le correspondía una dirección de máxima variación del campo escalar en el entorno 
del punto, representada por la del vector gradiente en esa posición de forma que dU = gradU : dí = V- dr. 


Vamos a ver ahora que para indicar la variación del vector de campo en el entorno del punto geométrico 
necesitamos “ascender” a todo un orden ¿Y = gFad V- dF = H: df por lo que recurriremos a la notación tensorial 
para ascender el orden del vector. 


ads 
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Como en el caso de los campos escalares y vectoriales, en general, supondremos que todas las funciones 
que intervienen en los campos tensoriales son continuas, derivables y uniformes, o lo que es lo mismo, le 
corresponde un sólo valor de cada función a cada punto del espacio, aunque la relación no será biunívoca pudiendo 
tener varios puntos un mismo valor para cada componente del tensor. Recurriendo, por tanto, al desarrollo en serie 
de la función vectorial AV = dV + > dV+..+ - d” Y y despreciando igualmente infinitésimos de orden superior 
tendremos: 


k 


¡+ 


7 LE cid p|_]9Y, ev, oy, A 9, o, => 9, or, A 
a? =dlv.i+v,) VAS 7 o e + el = dz|j + A rá = dz 


que agrupando en la forma: 


EAS AA A e 
Ox Ox Ox oy Oz lo) 
AAA BLAS laz 
ES 
Y El LAA CLASICA DE AS 2 loma 47 
Pa 7 > 2d [axi +dyj OE gradV- dr 


nos permite, en efecto, expresar la variación del vector del campo como “producto escalar” del “GRADIENTE DEL 
CAMPO VECTORIAL” en cada posición por el vector desplazamiento indicativo de cada dirección, o bien como un 
vector que tiene por componentes los productos escalares de los gradientes de las componentes del campo vectorial 
por el vector desplazamiento: 


> Vr Y >» 0 > 
—i+r—_j+—_k 
Oz 


; 2 Laxi+ dyj +daKJi+ 


Vr Vir Y > 
+ 2i+2j + —2k |: 
y? dz 


[axí + dy] +dz El. 


Y, vw, > wv > 
+ rs, kl: 


rado [dxí +*dyj +dek]k= 


= [rad V,: dr )í +|grad V,- dr |7 + [grad V,: dr | £ = grad V- dr 


y que no es otra cosa que el producto escalar de un tensor por el vector desplazamiento expresable más 
adecuadamente mediante la notación tensorial: 


Ox oy dz 
mn Y wo, dit E 
dV=|— — —]|- dy = gPad V : dr= HB dr (10.3) 
Ox dy Oz 
a 2 | a 
Ox oy Oz 


es decir, nuestro "vector" gradiente de un campo vectorial tendría 9 componentes a diferencia de las 3 del gradiente 
de un campo escalar, por lo que no se trataría de un vector sino de un tensor, como queríamos demostrar, que se 
denomina TENSOR DE ENTORNO ya que en una primera aproximación todos los fenómenos pueden considerarse 
lineales en este ámbito. 
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En este sentido, si consideramos el gradiente del propio vector de posición + obtendremos un tensor de 
particular relevancia por tener nulos todos sus términos salvo los de la diagonal principal:: 


ess 
dx oy 0 10.0 
> [Y Y Bl_ e 
VF = ss = 3” O 1 Ol=0 (10.4) 
Oz Oz Oz 0 0 1 
x y a 
denominado Tensor de Krónecker o tensor unidad. 
as OPERACIONES ENTES TENSORES 


Antes de seguir avanzando en el conocimiento de los tensores necesitamos previamente familiarizarnos 
con la estructura y la notación tensorial, para lo cual tenemos que definir las operaciones que admite este operador. 


- SUMA ALGEBRAICA 


Se entenderá que solo se podrá efectuar la suma algebraica entre tensores del mismo orden (cero, uno o 
dos) dado que, en una base determinada, cada componente del tensor suma será la suma de las 
componentes de igual subíndice de los tensores sumados: S; = T¡¡ + H, y se expresará como: 


Tales (10.5) 


Por el contrario, los productos entre tensores no imponen esa condición, pudiéndose efectuar entre tensores 
de diferente orden; es decir de orden cero (escalar), uno (vector) o dos (propiamente tensor). 


Distinguiremos entre dos tipos de productos ya que la diferencia fundamental entre ambos será que en 
unos (productos Internos) el orden del tensor resultante será la resta de los órdenes respectivos, mientras 
que en los otros (productos externos) el orden resultante será la suma de órdenes, para lo que resulta útil 
visualizar las rayas que indican el de cada tensor": 


PRODUCTO INTERNO O CONTRAÍD 


Este producto se obtiene por igualación de subíndices y suma entre los componentes de cada tensor, 
obteniéndose un tensor de orden resta el de los otros dos. Es el que sirvió de base para definir la 
transformada de un vector mediante un tensor, pero también lo es el propio producto escalar: 


W, 


VV=U 3 U=VYW 2D |Y,W]|:P,|= VW, + VW, + VW, (10.6) 


W, 


(m Es de notar que el producto de dos escalares (tensores de orden cero) es tanto interno como externo ya que (0 + 0=0-0=0) 


al igual que el de un es calar por un vector (1 + 0 = 1-0 = 1) y el de un escalar por un tensor (2 + 0=2-0= 2). 
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Ta Ta Tal ln | (7,1, +7,1,+T,V, 
Po V= W => VW. a TY > Dos Do y » V, = LY, + LY + LY, (0.7) 
Ty Ta Tal Y Ty Y, + T3V2 + Ta, 


T Ti, T € Ej) € 
1“ “1 * 5 
111€, +T19€12 +T,3€13 6 


tie=sE > E=tj, D ly Ty tl: [En En Es] = Tp *To7 HygEy + (10.8) 


+T¿¡€y, +T¿7 Ey, +T79€ 
3131 32732 “33-33 
Ta] Ty Ta € € €3 


Utilizándose en este último caso el punto doble punto por entenderse como una doble contracción. 


PRODUCTO EXTERNO O TENSORIAL 


Se define como un tensor de orden suma cuyas componentes son el producto de todas las componentes de 
uno por todas las del otro, que expresaremos en notación compacta por: 


oV=P > TaVe= Po (10.9) 


Sil 


VoHF=P > V,H, =P (10.10) 


siendo el primero un tensor de orden 3 que requeriría una representación en 3 dimensiones, mientras que 
el segundo, denominado PRODUCTO DIÁDICO de dos vectores resulta, efectivamente, un tensor. 


10.4 REPRESENTACIÓN GRÁFICA: CUÁDRICA TENSORIAL. 


A diferencia de los otros campos estudiados, no podemos encontrar una representación gráfica que permita 
visualizar globalmente la distribución de un campo tensorial por la dificultad que entraña el orden superior. 


Sin embargo, podríamos interpretar este orden 2 como un campo dentro de otro 
| campo ya que a cada posición en el espacio geométrico ordinario (variable 
y posicional F) se le asigna otro campo de variable direccional de forma que 

E estudiaremos primero cómo representar el tensor en un punto arbitrario para 


FER discutir posteriormente cómo se visualizaría el campo tensorial. 
PAY YY) y 
— AO YY) 


Supondremos, por tanto, el tensor T (referido a un sistema de ejes x;) asociado 
a un punto arbitrario que adoptaremos como origen de referencia dado que el 
Figura 10.2 vector transformado por el tensor tendrá el mismo punto de aplicación que el 


i=xyz 


vector original, 


Como se trata de encontrar una expresión analítica que identifique el tensor en cada punto, lo que estamos 
buscando es una relación entre el vector origen y el transformado o lo que es igual, una relación entre las 
coordenadas de los extremos de ambos vectores, por lo que elegiremos como vector origen el propio vector de 
posición + de ese punto P que representa el extremo de cualquier vector de forma que el transformado por el tensor 


> - A o ¡E 
será 7! - 7. 7 y Sus Componentes se podrán expresar como y! = To 
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Es de notar que el producto escalar de ambos vectores será un invariante ya que no dependerá del sistema 
de referencia elegido para su representación y por tanto, aunque cambien las componentes del tensor, como veremos 
en su momento, no lo hará cualquier magnitud escalar. 


En consecuencia, parece adecuado intentar representar espacialmente el tensor 
mediante uno de esos invariantes de forma que si buscamos el lugar geométrico 
de los puntos para los cuales el producto escalar de su vector de posición por 
el transformado mediante el tensor asociado al origen den un mismo valor nos 
encontraremos con la expresión: 


X, 


E y 
PF ap PS =x,T,x,= q (10.11) 


que desarrollada quedará: 
Figura 10.3 


Pu Zo TT [4% L 11% +7,9%, +7, 3% 


XX x0| Ta Ta Ta %|- 


2 2 2 
-= > + . 
q * * xa| Toy gp + gy | Ty yx + 227 HT ggx3 Ly +T ¡3% 30] +7 1X 7 +T 93037 +Ty1Xy%, +T39% 203 


La, Ta Tg|| Tax +Z 3x7 +T yx, 


y en el caso de tensores simétricos, para los que T¡¡ = T¡ (serán los que nos vamos a encontrar en nuestro estudio) 
la expresión obtenida se podrá simplificar en la forma: 


E 2 2 2 7 
q=T|¡X¡ + Tyx) + TygX3 +271,X,X, + 2T)3X,X, + 27, X3X, (10.12) 


ecuación de una cuádrica que, por tener como coeficientes los términos T; del tensor, recibe el nombre de 
CUÁDRICA TENSORIAL o forma cuadrática asociada al tensor y que corresponderá a un elipsoide, a un hiperboloide 
de una hoja o de dos hojas en función de los signos de los términos de la diagonal principal del tensor. 


Figura 10.4 


Como es de esperar, cuando los ejes de referencia coincidan con los ejes geométricos de las cuádricas, que 
denominaremos direcciones principales del tensor, la (10.12) podrá tener una expresión reducida en la forma: 


q =Ti¡2 + Ta? + Tx = 1 (10.13) 


al anularse los términos del tensor de subíndices diferentes y como el vector transformado será el mismo 
independientemente de los ejes de referencia adoptados, tanto el tensor asociado al punto como la cuádrica asociada 
al tensor serán únicos, variando exclusivamente su expresión. 


En consecuencia, generalizando a todo el espacio geométrico, la representación gráfica de un campo 
tensorial sería el conjunto de todas las cuádricas asociadas a cada punto del espacio, si bien solo tendrá utilidad la 
representación de algunas en particular, como veremos al estudiar el campo de tensores de inercia de un sistema 
material rígido, en especial las asociadas a los puntos de los ejes principales de la cuádrica correspondiente que, 
para nosotros, serán siempre elipsoides, denominados elipsoides de inercia del sistema material. 


De 
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Una de las ventajas inmediatas de la cuádrica tensorial se deduce fácilmente a 
partir de la expresión reducida (10.13) ya que el gradiente de la función q , normal 
a la cuádrica en cada punto de su superficie, como vimos en el tema anterior, 
tendrá precisamente la dirección del transformado del vector de posición del 
punto correspondiente: 


A Vg =2T,¡x,11, + 2T,2x,1l, + 2T,3x,81, = 2F' (10.14) 
Figura 10.5 


En efecto, si se han anulado los T¡, de diferentes subíndices, los ejes de referencia serán precisamente los 
ejes principales y el tensor habrá adoptado la forma diagonal por lo que el transformado de F será: 


pde 0 0 Xx 


ES A lo > po lx m0? 0. 7 o 
F"=x 4, +xU,+x34,=[ 0 T,j 0 2 | = T,x,4, + 77,x,8, + Tzgx3 4, (10.15) 


10) 
0 0 Tay 
lo que permitiría, mediante una construcción gráfica y el cálculo de su módulo según la expresión r'=-2_, 
r cosa 
determinar completamente cualquier transformado por el tensor, si bien su aplicación solo será práctica en el caso 
de tensores planos o mediante auxilio del ordenador. 


10.5 LIRECCIONES PRINCIPALES Y WALORES PROPIOS DEL  TENSOR, 


Como ya hemos intuido, si se eligen adecuadamente los ejes de referencia, la cuádrica asociada al tensor 
podrá tener una expresión reducida de forma que se anularán las componentes del tensor T;¡ de subíndices no 
repetidos i*+j quedando exclusivamente los elementos de la diagonal principal que habrán modificado su valor. 


La propia expresión del vector transformado por el tensor en su forma diagonalizada (10.15) nos permite 
dar una nueva definición de TENSOR como un “operador lineal que transforma cualquier vector multiplicando 
cada una de sus proyecciones sobre unas DIRECCIONES PRINCIPALES por un VALOR PROPIO de cada dirección”. 


Por tanto, a partir de esta definición, se deduce que cualquier vector situado sobre una de las Direcciones 
Principales del tensor, no modificará su dirección en la transformación, mientras que su módulo se verá afectado 
por el correspondiente Valor Propio o autovalor. 


Xa Denominando da, , d,, a, los unitarios según esas direcciones y A;, A), Ay los 
po. valores asociados respectivamente a cada dirección, la característica mencionada 
ed equivale a decir que el transformado de los propios unitarios se podrá expresar 


como: 
al =T-d=Ad 
Xo 1 1 1%1 
ms di=T:á,= Ad (10.16) 
2 2 NDA 
e A Ed 
dy = T* d, = Ad, 


de forma que la expresión del transformado de un vector referido a esa base 


— O0O> O> 10) A 
= + +v a, quedará: 
Figura 10.6 do: de Mat ds de To 9 


W=T-9=T-[v"a, +v/d,+v34)] =v Td, +v7T-d,+v3T:d,=v 4,4, +v, A,d, + vy A,d, (10.17) 


es 
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Si la base principal es ortogonal podemos expresar las componentes del vector referido a ella como 
o 


. . . (9) > > > > z . 
producto escalar del vector por el unitario correspondiente: v, =V: 4, =v; a, * d, y, análogamente, si la base 
arbitraria también lo es, los versores de las direcciones principales se podrán expresar en función de los cosenos 


> 


directores relativos a ella: a, = 0,1, de donde y,” = V:d, = v, 1, : 0, dl, = q,» Y sustituyendo en la (10.17): 


W o= A,v/ d, +A,v7 d, + A,v; 4, 

= A/A); U; +A7V¡Q7 07, ú, +A3V¡Qy ; 
= A,0Q,,Q,;V¡U; +A,0.7, 07, V¡U; YA¿ 03,03, V¡U = 
=[A, 0,0, *A,0),0) +A303,03] V¡U, = WU; 


Az 4, = 


(10.18) 


de donde se deduce que w, = [4, 0,0, +4, 0,,0, +4, Q,, 0] v, y al mismo tiempo w, = 7, v,, lo que nos 
permite expresar las componentes del tensor en una base cualquiera en función de los valores propios y los cosenos 
directores de los ejes principales del tensor con los ejes de la base elegida y la igualdad de los términos simétricos 
respecto a la diagonal principal: 


T;¡ =A,0,¡0,¡ *Az07Q7, ?A303/03 =T; (10.19) 


quedando condicionada la SIMETRÍA del TENSOR a la ortogonalidad de la base de referencia y de la Principal. 


Las 9 incógnitas de un tensor general se reducen, por tanto, en un tensor simétrico, a 6 , equivalente al 
conocimiento de los 3 valores propios, 2 cosenos directores de una de las direcciones principales y 1 coseno director 
de cualquier otra dirección ya que al ser ortogonal la base de referencia, la suma de los cuadrados de los cosenos 
directores del primer versor dará la unidad: % + d%, + o, =1 y al ser también ortogonal la base principal, el 
producto escalar de dos de ellas será nulo, proporcionando la ecuación 0, 0), + 0, ,, + 01, 0), = O, y la tercera 
se obtendrá como producto vectorial de las otras dos. 


10.6 DIAGONALIZACIÓN DEL TENSOR - INVARIANTES. 


Como ya hemos visto, el tensor es único para cada punto del espacio pero, al igual que los vectores, en 
función del sistema de referencia elegido así adoptará su expresión. Para obtener la forma diagonalizada, compuesta 
exclusivamente por los valores propios del tensor, es necesario encontrar, por tanto, la base principal y nos 
serviremos para ello de la propiedad definida en el punto anterior por la cual los transformados de los vectores 
contenidos en los ejes principales mantienen la dirección del vector original. En consecuencia, si wW = 17, el 
transformado de uno de los unitarios dá será: 


1] 
day lñ, 1, Lis] 2 %a = 7,0,+7,70,+7,30y Toda + Ty %+ 7 %=-=0 
A [Ta Ta Ta 2] > 1407 70) *T 070 y Ty + (3 )0+ Ta 0-0 (10.20) 
ns IZ, Lo Ta] | % May = 73/0,+73,0,+7330, Ty %* Ta %+* (y) =0 


sistema homogéneo de ecuaciones que para ser compatible tendrá que anularse el determinante de los coeficientes 


TA 7, Ta 
Tan Tak Ta |-0 (10.21) 
Za, Ta TyA 


obteniéndose en su desarrollo: 


¿Re 
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cri. 2.3 lu La Lo 
[AT ATA ata] [ne] [ln te Irrtil-o (10.22) 
11 22 33 T F T q T T 21 22 23 
y) 33 31 33 21 22 
Ta Ty Tx 


ecuación característica del Tensor cuyas soluciones corresponderán a los 3 valores propios que nos proporcionan 
la forma diagonalizada y que, para el caso de tensores simétricos serán 3 valores reales” 2, =4,, A,=4), Mk = 45. 


Ahora bien, dado que las soluciones de ) son unos valores escalares que no pueden depender de la base 
elegida para representar el tensor, se entiende que la ecuación característica proporcionará los mismos resultados 
independientemente de los valores de las componentes T; que aparezcan reflejados, por lo que los coeficientes de 
dicha ecuación tampoco variarán. Es decir, acabamos de obtener los 3 INVARIANTES más significativos del 
Tensor: 


1,= [2 +7 +73] = [4, +4, +43] TRAZA del TENSOR 
Ta Ta, T, T,3 Z, 1 2 0 1 0 | 1 0 
L= + E + + =4,4,+4A,+A,4, Y ADJUNTOS de la DIAGONAL PRAL. 
Ty, Ty La, Ty, a Ta 0 A, 0 A, 0 A, 
(10.23) 
TD To T,, lA, 0.0 
E, = Ta Tao Tas Al 0 A, 0 =A), A, A, DETERMINANTE del TENSOR 


T o 0.4, 


33 


pudiéndose haber deducido el primer invariante directamente de la (10.19) ya que tendríamos: 


2 2 2 
T,, = 4,0, +4,0, +Ay0%, 
2 2 2 
Ty = Aj02* A+ Ayo DT *Tyg* Ty = A(04,*077+ 043) =4,*A,*A, 


; 2 2 2 
Ty = A,0%3*A,0)3 +A30%7 


Finalmente, una vez obtenidos los tres valores propios del Tensor, las direcciones principales 
correspondientes se determinarían sustituyéndolos sucesivamente en las (10.20), si bien de las 3 ecuaciones obtenidas 
para cada valor de 4, una será combinación lineal de las otras dos debiendo completarse el sistema con la ecuación 
que proporciona la suma de los cuadrados de los cosenos directores de cada dirección, ya que partimos de una base 
ortogonal: dl, + dl, + o, = 1. Por eso mismo, obtenidas dos de ellas, podría obtenerse también la tercera recurriendo 
al producto vectorial. 


O En efecto, partiendo de una solución real 2 = A, (que siempre existirá) correspondiente a una de las direcciones principales á, 


tendremos que el transformado de ú, será otro vector 24, ; por tanto: 


UN Tú Tp Ti]! T,=>% »A»0 0 4 0 0 ak 
; 22 2 
0|= [Ta Tae Tafl9] > Tu=Tp=0 > [0 To To] > 0 Tagrk Ta |=(4-1 =0 
0 TT T 0 Tj= T,¿=0 Ts, Ty A 
a 32 3 ge e [0 7, Ty 0 Ty Tiyrhk 


y escribiendo de nuevo la ecuación característica: (T,, MT DT, T,=0 9 PAT) AT YAHT, Ti T4T,)=0 obteniendo las 
2 33 23” 32 a” 33 22” 33 23” 32 


soluciones: 


Ta TY TA TT, 7. 
hd EME ASE 2 Cato Tots) que para ser reales tendrán que cumplir (T,,+T,¿-4(T,¿Ty3"T,¿T,)20 


siendo, en efecto: Ti, *T5*2T Ty "4T y Ty +8T 93 T 93 TT] 2T 733 +8T yy T 37% (Ty "Ty +4T,¿T,, 20 ya que, en tensores simétricos, T,y=Ty 


y tienen el mismo signo además. 
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10.7 CAMBIO DE BASE. TENSORES ESPECIALES. 


Si lo que tenemos es una base ortogonal x”,x”,x”, variable (ejes móviles), y 
queremos referir los términos del tensor a dicha base a partir de los cosenos 
directores de sus unitarios respecto a otros iniciales (fijos) podremos escribir: 


1/= 0,1, * y 1,+ Uy, A tn [Ia 
ii¿= 0,8, + ay i,+0 8, > |1/|=|% % %llzl > 2'=43 (029 
ll, = 01, * Qyyti,* Oil, ll Az, Uy Ay [| a 

Figura 10.7 también expresables en forma reducida: ii, = 0%, 0 viceversa ¡, = 0,1; A 


— 


Por tanto, la relación entre las componentes de un vector referidas a ejes “fijos” Y = v,4, y a ejes 


. eS de e / E ES mn 6 / $ ce] 5 
“variables” Y = y, d, será v; =(Y:8,) = (v,11,* 0,,11,) = 0,, V, Y VICEversa: v, = 0, V;, O bien, utilizando la notación 


matricial, denominando A la matriz de cambio de base: 


3 _ = 
v/ 11 y %y Y (5) = 4 (5) 
/ a 
vw|=|%; %y Qallv] > y viceversa (10.25) 
/ 
Ya dj %z 03 [11 Y3 (9) =4 (51) 


expresiones igualmente válidas para los unitarios de las direcciones principales del tensor a,,d,,d, cuyas 
componentes designaremos ahora, para diferenciarlos de las de cambio de base como: 


de forma que las componentes del tensor referidas a los ejes fijos 7 ¡ =4,%,4, +4, 0), Ey +A3 0, €, =T) 
al referirlas a la base variable quedarán: 


* a) rd ad nn ar DI E > 4 ES Pa 
TAE 0 +A 0 +A 65 0 =T) T ¡5 TQ 0 
=A ¡Oy E yp0, E] HA7 0 oy 0 y HASO ¿AO 4 = =D y viceversa (10.26) 
E 4 7 7, 7. 0 = 5 dl % ¿Ol 
[A 4% +A7 Gp dy + Az Ep 831) 0,07 =T 900, Ti = Ty 0404 


= vil án* la ES 
y las componentes w, = Tv del vector vw transformado por el tensor, serán: w,* = 7, YT Oy Oz ly, Y 


A pesar de su complejidad la (10.26) puede expresarse en forma matricial si colocamos sus componentes 
en el orden correcto: 7, *= q,, 7, a, de forma que T"=ATA!: 
ij i l 


Ty Ty Ti Ay % %3|[7,, 7), Ti[[%, %, %, 
a To Ta |=[|% % %|lZ 72 Ta] [%, % 0% (10.27) 
T, 


a la Ta Az Ay Ay Ty, ES 


Az 03 Uy 


Los tensores que no cambian su expresión al cambiar la base se denominan INVARIANTES y constituyen 
un caso especial del que el más representativo es el tensor unitario o de Krónecker aunque, obviamente, todos los 
tensores de orden cero (escalares) también lo son. 
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IN PRINCIPIOS DE LA MECANICA 


En Cinemática hemos estudiado el movimiento como un problema aislado, sin tener en cuenta las causas 
que lo pudieran producir, pero el problema fundamental de la Mecánica consiste en encontrar o definir el estado 
de movimiento o equilibrio de los cuerpos conocidas aquéllas y para llegar al sólido rígido, idealización del sólido 
natural (deformable), es necesario previamente considerar una primera idealización al punto material. 


Si recordamos el primer principio de la Mecánica Clásica o PRINCIPIO DE INERCIA, "cuando sobre un 
punto no actúa ninguna fuerza, o está en reposo o se mueve con movimiento rectilíneo uniforme". Por tanto, 
en lo sucesivo consideraremos también éstas, es decir las FUERZAS, que definiremos como las causas que pueden 
originar el movimiento, modificarlo o incluso anularlo una vez establecido, como veremos más adelante y cuya 
relación expresa el 2 PRINCIPIO DE NEWTON o LEY FUNDAMENTAL, al postular la existencia de al menos unos ejes 
inerciales respecto de los cuales “la fuerza que actúa sobre un punto material le imprime una aceleración en 
la misma dirección y sentido que aquella” (fig.11.1), de tal manera que el cociente entre el módulo de la fuerza 
y el de la aceleración imprimida es una constante que se denomina masa inerte o inercial: 


F=ma > El=m (11.0 
siendo importante el cumplimiento de la condición de ejes inerciales porque de 
no ser así la aceleración absoluta, como hemos visto en Cinemática, sería 
d=d,+4d, y la (11.1) para un sistema de referencia que se trasladara con 
movimiento no uniforme, (además del término complementario en el caso de 


0 dicho movimiento no fuera de traslación) quedaría: 


Figura 11.1 F=md=m(d,+d,) = má, +mad, + F-má, =m 


de forma que la fuerza y la aceleración relativa a los ejes elegidos no cumplirían el principio, salvo que el término 
m d., se tradujese en una fuerza ficticia de inercia d = - má,. 


Esto nos lleva directamente al tercer principio, denominado PRINCIPIO de ACCIÓN Y REACCIÓN según el 
cual “cuando se aplica una fuerza aparece siempre otra igual y opuesta a ella que la equilibra”. 


O Si la fuerza se aplica a algo que no se puede mover la interpretación clásica será inmediata: aparecerán 
las fuerzas de reacción anulando las causas que lo tienden a mover: F+R =0. 


[a Si por el contrario nada impide el movimiento, aunque no vemos a primera vista esa fuerza que la 
equilibra, en algún sitio y de alguna forma tiene que aparecer, ya sea tomando ejes móviles ligados al 
cuerpo, lo que se traduciría en la fuerza de inercia Ú = - má de forma que respecto a ellos no hubiera 
movimiento al ser F + Ó = 0 (PRINCIPIO DE D' ALAMBERT), o bien, tomando ejes fijos de mayor "alcance" 
que incluirían el agente generador. 
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En este sentido el peso de los cuerpos, que siempre caemos en el error de considerarlo aislado, tiene sin 
embargo tiene su pareja en la propia Tierra o dicho de otro modo, en la fuerza con que el cuerpo atrae a la Tierra 
ya que no es otro que el principio de "atracción universal" (de acción de dos masas) que define los campos 
Newtonianos, según el cual: 


“las acciones que se ejercen entre dos cuerpos, reducibles a partículas 


E AR materiales (fig. 11.2), son dos fuerzas iguales y opuestas que actúan en la recta 
b, po CIZ, y Ha que les une, que son atractivas y cuyo módulo es directamente proporcional 
al producto de las masas! e inversamente proporcional al cuadrado de las 
y distancias”. 
E O [Fl =1F),1 = GuÉ (11.2) 
Figura 11.2 d 


n2 CAMPO Ss FUERZAS: CAMPO DE INTENSIDADES. 


Si nos limitamos a estudiar el efecto originado por una de las masas puntuales 
F=-G 2h  "m'p, definidas en el principio anterior vemos que, independientemente de la 
posición que pueda ocupar, la fuerza ejercida sobre m' es un vector dirigido 


m o r siempre hacia aquella, cuyo valor es función de la distancia que separa las dos. 
cd Esto responde obviamente a la definición de un campo vectorial de los 
estudiados en su momento ya que se podría expresar como F=F(Ff) y 

Figura 11.3 generalizando, a todo campo vectorial de la forma: 
F=F(F,1)=F(x,y,2,0) =F,(x,y,2, 1) P+F, (x,y 2,0] +F,(x,y,z, 0) £ (11.3) 


cuyos vectores representen fuerzas lo denominaremos CAMPO DE FUERZAS pero remarcando que se trata de un 
campo muy particular porque sería nulo en todos los puntos del espacio a excepción del que ocupase en cada 
instante la partícula material, 


Sin embargo dicho campo está "latente" y no solo se manifiesta cuando la partícula m' se coloca bajo su 
acción (fig.11.3), sino también para cualquier otra partícula m " que pudiera aparecer, ya que la fuerza con que m 
atraería »m " si ocupase la misma posición que m' sería un vector de igual dirección y sentido que la fuerza ejercida 
sobre ésta, cuya diferencia estaría tan solo en dicho valor. Podríamos decir por tanto que lo realmente existente es 
otro campo vectorial E = - G ml F,, asociado al punto geométrico y no al material y que la fuerza F que actúa en 
un instante dado sobre un punto material de masa m, representado por la posición concreta de la partícula (punto 
geométrico), será un vector F' = m E proporcional al vector correspondiente del campo F' que depende, por tanto, 


de los mismos parámetros que aquél. 


A este campo asociado al punto geométrico se le denomina CAMPO DE INTENSIDADES, aunque por 
cumplir dv = 0 se suele mal interpretar como el campo de fuerzas ejercido sobre la unidad de masa y a menudo, 
para designar la fuerza, en vez de F, se utiliza inadecuadamente f. 


(1 Se observará que en la (11.2) aparecen los términos y en vez de las masas m. Esto se debe a que, en primer lugar, queremos 
resaltar el concepto estático de este principio en contrapartida al concepto dinámico del 2” principio por el que la masa se definia 
implicando movimiento m = E y en segundo lugar porque aunque representan una medida de la materia, dimensionalmente no 
resulta muy coherente y obliga a dotar de unidades "residuales" a la constante G. 
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iz Supongamos una masa puntual situada en el punto y (X,,y,,Z,) (fig.11.4). En ese 
0) momento el vector de campo É ¡> Correspondiente a P,, representa la fuerza que 
pr SS actúa sobre ella. Si está en movimiento y describe una cierta trayectoria P, P, 
P, ..., la fuerza actuante irá variando, siendo en cada instante la fuerza del 
campo correspondiente al punto en que se encuentra el móvil, pero si permanece 

y en reposo seguirá actuando la misma fuerza del campo (si es estacionario) lo 
- cual parecer contradecirse con el segundo principio ya que no existiría 
aceleración, lo que indica la existencia de otra fuerza igual y opuesta con la 


Figura 11.4 misma recta de acción. 


La fuerza que actúa sobre el punto material tiene lugar, pues, por el hecho de encontrase éste en una cierta 
posición y poseer una masa sensible al campo de intensidades de que se trate. Es decir, la fuerza estará aplicada 
a la partícula pero su valor estará asociado a la posición respecto a la masa puntual central y existirá, 
independientemente de que esté o no ocupado por una partícula material. Si, por ejemplo, se trata un campo 
gravitatorio como el terrestre, éste producirá una fuerza (que denominaremos peso) atractiva sobre las masas pero 
no sobre las cargas eléctricas y en cambio, si se trata de un campo eléctrico, sucederá lo contrario, a no ser que se 
trate de un punto material sumergido a la vez en un campo gravitatorio y en uno eléctrico, sufriendo los efectos de 
ambos sometido a la resultante de las dos fuerzas activas. 


La más importante de las fuerzas que conocemos, el peso, no es otra cosa que la fuerza con que la Tierra 
atrae los cuerpos que se encuentran bajo los efectos de su campo de intensidades F' = mE. De ahí que "la 
gravedad" no sea ni una fuerza ni una aceleración, como se dice erróneamente muy a menudo, sino una intensidad 
y como tal sus unidades son [F/m] a pesar de la similitud de dicha expresión con la del 2* Principio: F' = má. 


La diferencia entre E y d es que la primera es la causa, la intensidad existente en esa localización 
independientemente de que haya masa o no, mientras que la segunda es el efecto, lo que le pasa a la masa que se 
introduce en el campo y de no existir campo (causa) no le pasaría nada (efecto) a la masa. Es, en consecuencia un 
ejemplo muy ilustrativo de vectores idénticos ya que £ = d y tienen el mismo punto de aplicación. Sin embargo el 
significado de ambos vectores es opuesto: causa = efecto o, dicho de otro modo: la intensidad del campo en un 
punto sería igual a la aceleración que imprimiría a una masa que estuviera localizada en él. 


m3 CAMPOS CENTRALES. CAMPOS CONSERVATIVOS 


En cualquier caso, la definición dada para los campos denominados NEWTONIANOS representa un caso 
particular incluido en otro grupo mucho más general que abarca todos los CAMPOS CENTRALES, cuyos vectores 
tienen por dirección la que une el punto de aplicación M con otro punto fijo O que se denomina centro del campo: F' =1F 
y por módulo una función de la distancia entre ambos || = F(»). 


De este modo, si + es el vector de posición de la partícula material adoptado el 


le ] 
| F punto fijo como origen de coordenadas, (fg.11.5), llamando F, al vector unitario 

l— qe > en esa dirección podremos escribir el vector campo en la forma general: 

o P 
ZO y = F(r) . UR, (11.4) 
v Como otro caso particular de los centrales, por tanto, los campos Newtonianos 
Figura 11.5 de la (11.2) serían de la forma |F| = F(r) = Z. 
di 
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Pero vamos a encontrar otra aplicación inmediata en las fuerzas originadas por 
las deformaciones lineales elásticas cuyo comportamiento se asemeja al de un 
muelle  (fg.11.6), al ser el módulo proporcional al alargamiento 
|F| =k Ar =k(r - 5) = F(r), donde 5 representa la longitud natural, tener la 
dirección del muelle y sentido hacia el extremo del muelle o el opuesto según 
sea el de recuperación: F'= - kArF,. 


Figura 11.6 


Por último, dada la importancia de los campos centrales y la asiduidad con que nos los vamos a encontrar 
en Mecánica, recordaremos una de sus propiedades más importantes: la de que todo campo vectorial central es 
CONSERVATIVO, es decir, deriva de una función escalar de forma que F=E(r) F, = grád U. Una forma de 
averiguar si un campo vectorial es conservativo o no era comprobar la igualdad de las derivadas cruzadas de sus 
componentes. Pero había otra forma mucho más inmediata que era encontrar la función U de donde procedía ya 
que, de existir, 

F"=gradU > dU=grádU:dr=F"-dr => U= f[F":dr (11.5) 


y puesto que todo campo central es función de una sola variable, r, resulta inmediata dicha integración: 


U= JP de= [Poy de [Fo d PLr k 


Fm ar]=[Fom2 ajr?] = [ Foyar = 7% 
[OF ars [Far = [Far =U(») 

En consecuencia, las superficies equiescalares de los campos centrales serán esferas. 
tia TRABAJO. ENERGÍA. POTENCIAL 


En principio la idea de trabajo implica "algo" que lo realiza. Supongamos pues un campo de fuerzas, nulo 
en todos los puntos menos en aquellos donde se localice materia, y una masa puntual m en dicho campo que, por 
efecto de la fuerza a la que va a estar sometido adquiere una aceleración y una velocidad que origina un cambio 
de posición de un punto a otro del campo donde otra fuerza le imprime otra aceleración y, en consecuencia, 
modifica la dirección de la trayectoria ... ; es decir, estamos considerando una trayectoria real. 


El primer punto al que tendería a irse la masa de forma natural por efecto de la fuerza actuante del campo 
sería uno de la trayectoria real, realizando, por tanto, un DESPLAZAMIENTO REAL dF NATURAL. Pero podemos 
tambien "imaginar" un movimiento cualquiera de la masa con independencia de la fuerza actuante que no tiene 
por qué ser el que realizaría de forma natural y lo llamaríamos, por tanto, DESPLAZAMIENTO VIRTUAL 07. 
Imaginemos, por último, la masa puntual "ensartada” en un alambre; el único desplazamiento que "podría" realizar 
de entre los que pudiéramos imaginar sería a lo largo del alambre, es decir el único REAL COMPATIBLE CON LOS 
ENLACES, independientemente del que seguiría libremente bajo la acción del campo de no verse forzado en su 
trayectoria. 


Bajo estas premisas y recordando la definición de circulación de un vector de 
campo a lo largo de una línea, definiremos el trabajo elemental que ejercerá el 
punto material para pasar de un punto geométrico a otro (lo suficientemente 
próximo como para aceptar que el vector del campo de fuerzas apenas se ha 
modificado) como el producto escalar 58 = F- 8 si el desplazamiento es 
virtual 6 dE = F- dr si el desplazamiento es real, denominando al primero 
trabajo virtual y al segundo trabajo real, siendo éste uno sólo de todos los 
virtuales posibles. 


Figura 11.7 
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El concepto de "imaginar" por donde se vaya a romper una cosa es importante bajo el punto de vista de 
los trabajos virtuales, ya que cubre todas las posibilidades, pero su aplicación está limitada al estudio en el entorno 
de los puntos ya que no existe trayectoria. 


Sin embargo el Trabajo total de la partícula para pasar de un punto 1 a otro 2 (£g.11.7) sólo tiene sentido 
como suma de todos los trabajos elementales reales a lo largo de los puntos de la trayectoria real, que no tiene por 
qué ser la natural (la que originaría la fuerza) sino venir forzada por algún tipo de ligadura: 


2 
E = [Fd (11.7) 
1 


Esta expresión se parece a la de la circulación de un vector, pero solo formalmente, no conceptualmente, 
ya que la circulación es un concepto abstracto que, al no implicar materia, no puede conllevar qué cosa es la que 
realiza trabajo. Solo si existe materia podemos hablar de fuerza y, por tanto, de trabajo. Pero, relacionándolo con 
la circulación, podemos intuir que si la partícula pasa de 1 a 2 por otro camino no tiene, en principio, por qué 
realizar el mismo trabajo. Sólo ocurría esto en el caso de que el campo de vectores fuera conservativo (procede de 
la aplicación del gradiente a una determinada función escalar). 


En tal caso, si el campo de fuerzas es conservativo, diremos que procede de una cierta función escalar 
U denominada función de fuerzas que, por conveniencia, definiremos como la función de Energía Potencial 
cambiada de signo: £,, = - U de forma que: 


E = grád U = - grád E, (11.8) 


y el trabajo resultará: 
S á - cd M_ 0) 
8 = rádU + de = [ -grádE. + dr = - dE. = Bx*=E 
; did 5-[ g / RE / (1.9 
E 


no dependiendo del camino seguido, sino de los valores que asume en los extremos de la trayectoria la función 
E, =- / F (P) di = E, (*) denominada de ENERGÍA POTENCIAL, es decir, el trabajo realizado por la partícula 
sometida a las fuerzas del campo conservativo es igual a la variación de energía potencial entre el origen y el 
punto final, de forma que el trabajo es positivo cuando se deja llevar de forma "natural" por ellas ya que los puntos 
tienden a ir de mayor a menor energía potencial. O dicho de otro modo, si la fuerza que hace moverse al punto (no 


siempre la que actúa) es la que realiza el trabajo, éste será positivo por tener la dirección del movimiento. 
Procediendo de forma análoga con el campo de intensidades E = £ E. obtenemos la equivalente a la (11.9): 


2 2 2 2 Y, 
Le? | = | Esdi= [gradl-di= | -grád od - [ogridv-d=- [| dV=V,-V, a 
: 1 1 1 Y, 

en la que se ha denominado POTENCIAL o Función Potencial a la energía potencial por unidad de masa 

V= = E - fÉ E* dr = V (F) de forma que, aun no teniendo sentido hablar de trabajo por no tratarse de una fuerza, 

la (11. 10) se podría entender como el trabajo realizado por la masa unidad al cambiar de posición, cuyo valor vendrá 

dado por la diferencia de potencial entre la posición inicial y la final, de forma que si ambas se encuentran en la 

misma superficie equiescalar (ahora denominada EQUIPOTENCIAL), no se realiza trabajo independientemente del 
camino seguido. 


Ahora bien, si queremos averiguar en qué se invierte la energía en los campos conservativos tendremos 
que introducir la relación causa (F) > efecto (m d) a nivel de trabajo, para relacionarlo con el movimiento de la 
partícula material. 
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ms TEOREMAS FUNDAMENTALES DE LA MECÁNICA 


TEOREMA DE LA ENERGÍA Cinenca | 


A partir de la definición de trabajo de una fuerza y conectando con lo estudiado en cinemática podremos 
expresar el trabajo elemental de las fuerzas aplicadas sobre un punto material en la forma: 


EA PA A 2 [mv 
di di ai 2 


=<l 
x 


= dE, (1.19) 


donde se ha denominado energía cinética al término . jp y?. 
Z 


Esta expresión del Teorema de la Energía Cinética, por la cual el trabajo elemental realizado por las 
fuerzas equivale a la variación de la energía cinética del punto material, es válido para todo tipo de fuerza 
aplicadas a la partícula y constituye una magnitud escalable por lo que, de estar aplicadas varias fuerzas, la 
variación de energía total sería la suma de todas las producidas por el trabajo elemental de cada una de las fuerzas, 
cuyo signo podría ser negativo en función del ángulo que formase la fuerza con el desplazamiento elemental. 


Es, en definitiva, la ecuación analítica que va a regir la Mecánica, en términos de trabajo, como 
alternativa al planteamiento en términos de fuerzas, basado en la ecuación vectorial. 


PRINCIPIO DE CONSERVACIÓN DE LA ENERGÍA 4 


Ahora bien, si TODAS las fuerzas aplicadas a la partícula derivasen de un campo escalar se cumpliría, por 
lo visto anteriormente, que dé = - dE, = dE, y el trabajo total realizado por las fuerzas del campo para pasar 
de un punto a otro sería: 


2 
1 2) 0 a 
4 = EV -E9 = EP - EY (11.12) 


es decir, podríamos afirmar que si el trabajo es positivo lo que se gana en ENERGÍA CINÉTICA (efecto) se pierde 
de ENERGÍA POTENCIAL(causa) y viceversa. Por eso se utiliza la función de Energía Potencial E, en lugar de la 
función U = -E, ya que en tal caso tendríamos que decir que lo que se gana en una se gana en la otra y no da la idea 
de inversión. 

Asociando términos en la (11.12) para relacionar lo que pasa en la posición 
inicial con lo que pasa en la final (Ag. 11.8) nos queda: 


0 pa AS e ES +ES (11.13) 


expresión del PRINCIPIO DE LA CONSERVACIÓN DE LA ENERGÍA MECÁNICA 
por el cual, se mueva el punto como se mueva, si la fuerza que lo hace moverse 


Figura 11.8 
deriva de una función escalar, se conservará la energía mecánica total. 


De ahí que se denominen conservativos estos campos, en los que se verifica, como se recordará, 
rót F?* =0, en contraposición a los disipativos cuyo rót F'+ 0 y en los que la energía mecánica no se conserva sino 


que tiende a disminuir, pasando a energía de otra naturaleza, muy frecuentemente térmica. 
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Dada su gran importancia para el desarrollo de los teoremas que permitirán abordar el estudio de los 
sistemas materiales y, en especial el del sólido rígido, veremos también los teoremas relativos a la CANTIDAD DE 
MOVIMIENTO, definida como magnitud vectorial: m Y que constituye la clave de los métodos alternativos basados 
en la expresión analítica de la energía mecánica del punto material. 


TEOREMA DELA CANTIDAD DE MOVIMIENTO | 


A partir de la ecuación fundamental definida por el 2* Principio podemos hacer 


A + dv _ d(mv 
Rsni=n Ll . PRD (11.14) 
dt dt 
dd a relacionando directamente la fuerza aplicada al punto material con la derivada 
Figura 11.9 respecto al tiempo del vector CANTIDAD DE MOVIMIENTO que acabamos de definir. 


La expresión anterior, en la forma ¿7 = Fdt = d(mv), nos define el IMPULSO elemental de una fuerza 
entendida como la variación de la cantidad de movimiento que imprime en el punto material. Por tanto su 
integración, que no tiene por qué ser fácil si tenemos en cuenta que la fuerza puede ser F'= F (y, F' £), nos dará el 
impulso total de la fuerza como incremento de la cantidad de movimiento entre la posición inicial y final 


F = [Far = Amv (11.15) 


que, supuesta a partir del reposo, permite diferenciar claramente entre fuerza e impulso de una fuerza al ser la 
primera un vector proporcional a la velocidad alcanzada 7 =my y el segundo un vector proporcional a la 
aceleración: F' = má. Casos particulares serán: 


F =cte => IT = FAt= AmV = cle 
P=), > 1 =FAt=Amv (mov. rectilíneo) 
F=0 > I1=0 => mv =cle 


TEOREMA DEL MOMENTO  CINÉTICO 


Por definición de momento de un vector, diremos que el MOMENTO de una 


_ FUERZA respecto a un punto A (fijo o móvil) será el producto vectorial 
mv 


> — dv d(mv) 
M,=BxXF =B6xmd = 6xm— = $x —_—< 
¿=D pxm pxm Pr p Fr (11.16) 


Va mientras que denominaremos MOMENTO CINÉTICO al momento de la cantidad 


“o  ""— — demovimiento respecto a dicho punto: 
9, = PxmvV (11.17) 
Figura 11.10 
Por tanto, si derivamos esta última expresión, teniendo en cuenta que si A no es fijo tendrá velocidad 
; PEÓN dp de dr, 
respecto a un punto fijo O de forma que fp =F-F, y e ide e _ = Y — V, tendremos: 
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dó P ? 
A dmv = (9-7) xmv+ pr 


A = M,¿-Y,xm (11.18) 


= Depor 
É 


Solamente si A es un punto fijo, por ejemplo el origen de referencia O, o si se mueve en todo instante con velocidad 
paralela a la del punto material el término v,x mV =0 y la expresión del momento cinético quedará: 


M,=-— (11.19) 


por la cual “la variación del momento cinético de un punto material respecto a un punto fijo o cuya velocidad 
sea en todo instante paralela a la del punto material es igual al momento, respecto a dicho punto, de las fuerzas 
actuantes (activas y pasivas) sobre la partícula material”. 


Independientemente de la importancia que van a alcanzar estos teoremas en la Mecánica del sólido rígido 
es de remarcar el paralelismo entre la (11.14) y la (11.19) ya que la primera establece la igualdad entre la fuerza y 
la derivada del vector cantidad de movimiento, mientras que la segunda lo hace entre el momento de la fuerza y 
la derivada del momento del dicho vector. 


nó ECUACIONES DEL MOVIMIENTO PARA EL PUNTO MATERIAL 


A la vista de los Teoremas analizados podemos, cuanto menos, asegurar que no siempre habrá un “único” 
planteamiento posible para resolver los problemas de la Mecánica, pudiéndose diferenciar claramente dos 
procedimientos alternativos de estudio, según se adopten planteamientos de trabajo-energía (escalares) o de fuerzas 
(vectoriales): 


ENERGÍA dE = F «di = dE, ysi F=-grádE, > d£ = dE, = - dE, 
dudis Pa ó sa a» (1.20) 
FUERZAS P=-YEF =máa=mts. te o, si el punto es ligado,  M,= As po 
di di o «e di di 


si bien los valores que puedan determinarse por ambos, obviamente, no van a diferir. 


El planteamiento dependerá fundamentalmente del objetivo ya que por el primero se podrán determinar 
valores escalares, y recurriremos a él si deseamos conocer solamente las condiciones particulares en una 
determinada posición, incluida la posición de equilibrio, pero no permiten por sí solas tener definido el movimiento 
en su totalidad salvo en los casos en que la trayectoria del punto venga limitada por algún tipo de ligadura.. 


Como fin último, el objetivo de la Mecánica es poder conocer en todo instante la posición de los cuerpos 
sometidos a un sistema de fuerzas determinado y la del punto material, según vimos en Cinemática, podía quedar 
determinada indistintamente mediante el conocimientos de alguno de los siguientes datos: 

O lasecuaciones paramétricas del movimiento F=F(t) 
1  latrayectoria I' y la ley de arcos s =S (t) 


1D  latrayectoria I' y la hodógrafa E ... 


Por tanto, en una primera aproximación distinguiremos dos tipos de problemas según se tengan definidas 
inicialmente las causas del movimiento o el movimiento en sí que, por otro lado, no es lo que habitualmente se va 
a conocer. Es decir, el problema se resolverá: 
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o 


1 estudiando el movimiento de la partícula F =F(t) que al cambiar de posición va sufriendo las fuerzas 
E 2, Ñ 
F=F(?, t) correspondientes a la aceleración á(F,t) = = en la posición que ocupa en cada instante 
a di 
má(*?,t) = F(?,t) (problema inverso). 


Ln estudiando las fuerzas perseguidoras del punto material que modifican la velocidad con que lo mueven 
Le F determinando su posición en cada instante en función de unos datos iniciales que determinan las 
m 
constantes de integración ña = E (problema directo) 
di ”m 


En este último caso, conocidas las fuerzas, la integración sucesiva de las ecuaciones que proporciona el 
problema, con las condiciones iniciales de posición y velocidad, nos lleva a la obtención de las ecuaciones 
paramétricas del movimiento que, en cartesianas expresaremos como: 


dl A Xx =Xx(0) +v,, x= x(1) +v,t+x, 
j=-=F, [2 I=30+v), [3 YI =I0 +vgtty, (11.21) 
Pel 2 = 2(t) +v,, 27 =2(() +v,t+z, 


aunque en ocasiones, como en los campos centrales, pueda sernos más útil expresar el movimiento mediante 


coordenadas polares F = rF, ya que el plano osculador de la trayectoria será el mismo en todo instante al tener 
la aceleración la misma dirección que el vector de posición y, por tanto, el movimiento será plano: 


A > v =P > r=r(b) 
á dro y m á 


: » . ; (11.22) 
a)=2+0+r0=22(0)=0 >  v=r0 > 0=0( 


o bien, cuando se conozca la trayectoria, refiriéndolo al sistema de ejes intrínseco al movimiento, es decir, el 
formado por la tangente, la normal y la binormal: 


1 _ dy 
o v=v(t == s=s(f) 
n= lr, = 2 fa fo (11.23) 
ag = 2F,=0 p=p0 == Y = Pin 


En cualquier caso, resulta imprescindible para su resolución el conocimiento de los valores que pueden 
adoptar las fuerzas en el trascurso del tiempo y su posible dependencia de otros parámetros (función de la posición 
del punto, del tiempo, de la velocidad ... etc. F =F(v,F, f)) siendo de particular importancia para nuestro estudio 
las fuerzas que dependen exclusivamente de la posición F'=F(F) y no de las otras variables como el tiempo 
(campos estacionarios) o la velocidad. 


17 EQUILIBRIO DEL PUNTO LIBRE. ESTABILIDAD 


Entendiendo como estado de equilibrio del punto material aquél en el que se den las condiciones 
necesarias para que no cambie de posición en el transcurso del tiempo respecto de un marco inercial, o lo que 
es lo mismo, que no tenga velocidad v=0 y, lo más importante, que tampoco la tenga en un instante anterior o 
posterior, a=0, diremos que la ESTÁTICA será la parte de la Mecánica que estudiará las circunstancias para que 
pueda darse el equilibrio así como las características del mismo, mientras que la DINÁMICA se entenderá como la 
parte de la Mecánica que determinará las características del movimiento. 
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Por tanto, como caso particular de la Mecánica, ambas seguirán planteamientos paralelos, entendiendo 
que la partícula estará en equilibrio en una posición determinada, en términos de fuerzas, si sobre ella no actúan 
fuerzas o se anulan entre sí: 3'F"=md'=0 (ya que en tal caso a =0 <> dv =0) 0 si el momento de las fuerzas 
respecto a un punto fijo es nulo (ya que entonces mv'=c” 2 dv =0) mientras que, en términos de trabajo, 
tendrá que anularse cualquier trabajo elemental de las fuerzas aplicadas a la partícula en el entorno de la posición 
de equilibrio (trabajo virtual, no real) 4 £=3F, :9F"=0 (ya que d £ =m v dv =0) y, por lo mismo, si las fuerzas 
derivan exclusivamente de un potencial, que éste alcance valores extremos (máximos o mínimos) dado que 
dé =-dV =0. 


En ambos casos se ha impuesto la condición dv =0 ya que si a=0 también se cumplirá a, =0 y 
dv = a, dí = 0 pero también debe presuponerse que la velocidad en cualquier instante es nula ya que podría estar 
moviéndose con velocidad Y constante (movimiento rectilíneo y uniforme) y carecer por tanto de aceleración desde 
el punto de vista del planteamiento de fuerzas, o análogamente con movimiento circular uniforme y tener 
componente normal la aceleración pero cumplirse que el trabajo es nulo por serlo dv pero no y. 


En definitiva las ecuaciones del equilibrio de un punto material libre se pueden resumir en el siguiente 
esquema que se reproducirá en lo sucesivo tanto para el punto ligado como para los sistemas materiales y en 
particular el sólido rígido: 


EF, =R=má=0 Fuerza resultante 
e ds E. Clásica 
7 e =PX 7 =0 Momento resultante 
(11.24) 
58 = EF 5P=0 Trabajos virtuales 
E. Analítica 
dé =-YdV =0 Potenciales 


ESTÁTICA CLÁSICA 


En definitiva, si la condición de equilibrio es que el punto no cambie su posición a lo largo del tiempo 
respecto de un marco inercial, adoptando un sistema cartesiano de referencia el problema se resolverá 
determinando el punto para el que F(*) = 0 aunque en otros puntos de su entorno no se cumpla (sin 
olvidar que la fuerza actuante puede tratarse de la suma de varias fuerzas F ¿) y como no puede depender 
del tiempo porque el equilibrio no debe romperse con su paso tendremos 
F(xyz)=0 = F(xyz) ia Fxyz)j +F (xyz) k, lo que nos proporciona un sistema de tres ecuaciones 
con tres incógnitas (x y z) coincidentes con los grados de libertad del punto libre en el espacio, cuya 
resolución nos determinará, de existir, la posición de equilibrio (aunque encontraremos más de una en 
general): 

F(xyz)=0 X=X, 

F(xy2)=0 3 y=y, (11.25) 

F (xyz)=0 z=2Z, 


La ecuación de momentos resulta trivial en el caso del punto libre puesto que si la resultante es cero y el 
momento en un punto también lo es, al considerar el que ocupa la propia partícula donde están aplicadas 
las fuerzas nos remite a la misma ecuación de las fuerzas. Sin embargo, su utilidad comenzará a ser 
patente e incluso mayor que la de fuerzas en el punto ligado. 
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ESTÁTICA ANALÍTICA ] 


«ds 


TRABAJOS VIRTUALES: En esencia se trata de encontrar la posición para la cual las fuerzas 
actuantes sobre el punto van a equilibrar sus trabajos teniendo en cuenta que unas darán trabajo positivo 
y otras negativo para un desplazamiento virtual. d 4 = F, :9r'+ F,:9F'+ F, :ÓF... =0 


Como 5 se puede descomponer, para cada fuerza considerada, en una 
componente según la dirección de la fuerza y otra normal, que no daría 
trabajo por tratarse de un producto escalar, basta con plantear 
desplazamientos en la dirección de cada una de las fuerzas y sentido de 
incrementar las variables correspondientes para simplificar el 
problema, teniendo en cuenta que en la expresión final del trabajo 
elemental no podrán aparecer más variables independientes que grados 
de libertad. 


Por ejemplo, supuesto un punto pesado sometido además a 3 fuerzas 
centrales que dimanan de los puntos O, A y B, (la primera atractiva y 


Figura 11.11 las otras dos repulsivas) directamente proporcionales a la distancia y 
de igual coeficiente k, haríamos: 
2 
ro= x2+ y? r.Sr,= x8x + y0y SS =- kr Sr, + kr Sr, + kr,Or, - PSy = 


rSry==(a-aJóxtyóy > =-k[x8x+y8y]-k[(a-x)6x-y8y]+ k[x8x-(b-y)8y]-P3y = 
rjOry= x0x - (b-y)0y = k[x-a]8x + [k(y-b) -P]8y = 0 


ri=(a-xpry? > 


rg =x2(b-y 


de forma que para anular el trabajo tendrán que anularse los coeficientes de cada una de las variables x 
e y ya que el punto tiene dos grados de libertad, obteniendo la posición de equilibrio en x=a, y =- +b. 


POTENCIALES: Un caso particular notable lo van a proporcionar los campos conservativos 
F* = - grád E, dado que: 


* __9E, _ 
F, id 
» __9E 
E A DS E,(xyz)=E, (11.26) 
F¡=-Lz=0 
dz 


es decir, la POSICIÓN DE EQUILIBRIO de un punto sometido a un campo de fuerzas conservativo 
coincide con los extremos de la función de energía potencial y en tal caso solo necesitamos conocer la 
función E, para obtener la posición de equilibrio sin más que anular sus derivadas parciales respecto a las 
variables independientes x, y, z o su diferencial exacto en caso de que dependa de una única variable (p.ej. 
los campos centrales) lo que equivale a plantear el equilibrio en términos de trabajo ya que: 


dE (0) _ 
dr 


|gradE, | > dE,(r) = grád E, df =-d8 =0 (11.27) 


En el ejemplo anterior, suma de cuatro campos conservativos, la expresión quedaría exactamente igual: 
dE,> kr,¿Sr, - kr Sr, - kr Sr y + POy = k[x8x+y8y] + k[(a-x)0x-y8y]-k[x8x-(b-y)8y]+P8y = -k[x-a]óx - [k(»-b) -P]8y = 0 


ya que, obviamente no es necesario encontrar la función de energía para volver a derivar. 
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Como los valores extremos de una función pueden representar los máximos, los mínimos, los casi máximos 
y los casi mínimos y veíamos que un punto tiende a ir naturalmente de donde hay mayor energía potencial a donde 
es menor, podemos deducir la importancia de este factor en la configuración del equilibrio, que nos va a permitir 
asimismo determinar la Estabilidad. 


El concepto de ESTABILIDAD DEL EQUILIBRIO es tanto o más importante que el equilibrio mismo 
ya que, en realidad, las cosas no es que vayan a estar siempre quietas, sino que van a moverse poquito. Por ejemplo, 
una gradería vibrará por acción de las cargas móviles, pero al ser estable lo va a hacer imperceptiblemente, 
pareciendo casi quieta. De hecho, las Estructuras se denominaban tradicionalmente Estabilidad de las 
Construcciones. 


En ese sentido diremos que el Equilibrio es ESTABLE en un punto cuando al desplazarse ligeramente (en 
cualquier dirección) de la posición de equilibrio (4'(F,) = 0) las fuerzas del campo correspondientes a la nueva 
posición (F(F, + 8F) + 0) (un poco diferentes por lo general a no ser que se trate de un campo uniforme o 
descontinuo) tiendan a devolverlo a su posición inicial. 


Análogamente, el Equilibrio en un punto será INESTABLE cuando al desplazarse ligeramente de la 
posición de equilibrio las fuerzas del campo correspondientes a la nueva posición tiendan a alejarlo aún más de su 
posición inicial y, por último, será INDIFERENTE cuando al desplazarse ligeramente a cualquier otra posición del 
entorno de la posición de equilibrio, ésta siga siendo de equilibrio. 


Dado que la posición de equilibrio corresponde a los valores extremos de la 
función de Energía Potencial, al desplazarnos a cualquier otro punto de su 
entorno, si se trata de un mínimo será E,>E,(e) y como el punto tenderá a ir de 
donde es mayor la energía potencial a donde es menor, volverá a la posición 
inicial luego es estable. Si está en un máximo: E,<E,(e) y se alejará más aún 
buscando otros puntos de menor energía potencial mientras que si E,=E,(e) el 
Figura 14.12 punto seguirá en equilibrio. 


En general es muy complicado plantear el estudio de la estabilidad en campos no conservativos ya que 
implica analizar los valores que adoptan las fuerzas para pequeños desplazamientos en cualquier dirección del 
entorno de la posición de equilibrio. 


Sin embargo para los casos particulares de campos conservativos es inmediato ya que tendremos 
únicamente que analizar la variación de la función escalar E, y determinar si aumenta o disminuye su valor, lo que 
se consigue PE el pl de las derivadas sucesivas ya que por el desarrollo de Mc Laurin: 
E, 7 Ey(e) = dE, + > PE +, d”E, Cuando hay equilibrio el primer término es sistemáticamente nulo y el 
signo de la Aieenda lo da el término de la derivada segunda o, como caso particular en que también sea nulo, se 
verá la tercera. 


Mediante la exposición anterior hemos obtenido el denominado TEOREMA DE LEJEUNE DIRICHLET, que 
establece que, en las posiciones de equilibrio, el potencial pasa por un valor extremo; siendo inestable si dicho valor 
es un máximo y estable si es un mínimo. 


Supongamos que uno de estos puntos, P, corresponde a un máximo. En tal caso, para los puntos próximos 
aP el gradiente de V estará dirigido hacia ese punto y, por tanto, É = - gñad V en sentido opuesto. Esto significa 
que, situado el punto material en P, se encontrará en equilibrio pero, al desplazarlo ligeramente de esta posición, 
actuarán sobre el punto fuerzas cuya resultante tenderá a separarlo aún más, tratándose, en definitiva, de un 
equilibrio inestable. 
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Si, por el contrario, V fuera un mínimo en el punto P, cualquier pequeño desplazamiento de la posición 
de equilibrio daría lugar a fuerzas cuya resultante tendería a llevar de nuevo el punto a dicha posición, lo que define 
un equilibrio estable. 


Si, en una región del espacio, el potencial V se mantiene constante, en ella se verificará R=0 y será 
posible desplazar el punto de un lugar a otro de esta región sin que aparezcan fuerzas que favorezcan o se opongan 
a su movimiento, En tal caso el equilibrio se denomina indiferente. 


Si la función depende de una sola variable su determinación será inmediata. En caso contrario se tendrá 


que recurrir al Hessiano. Si Y =V(p,n) habrá de cumplirse > = 5 =0) y para determinar la naturaleza del equilibrio 
estudiaremos el signo del hessiano: 
PY 3Y 
op? pon 
H = 
PY er 
onóp an 
de forma que si: 
2 
Lo máximo 
mp. * 
2 
Y mínimo 
op? 


H<0 no hay máximo ni mínimo 
H=0 - debe proseguirse la investigación. 


Análogamente se procedería si tuviéramos el caso en que Y =V(p,n.v). 


18 PUNTO LIGADO. EQUIVALENTES DE LAS LIGADURAS 


Obviamente las fuerzas (las de los campos) serán las causas que podrán producir el movimiento del punto 
material pero la experiencia demuestra que, aun existiendo fuerzas actuando sobre la partícula, cuando se da algún 
tipo de ligadura el movimiento no siempre se va a producir. por lo que las ecuaciones del movimiento no 
responderían a la realidad. Necesitamos, por tanto, definir unos equivalentes que “revaliden” las ecuaciones del 
movimiento o la energía para el punto ligado, es decir, cuando tenga restringido alguno de sus grados de libertad 
que, en el espacio ordinario, coincidirán en número con las coordenadas necesarias para determinar su posición. 


Esta limitación podría venir originada al obligarlo a moverse por una superficie 
determinada, quedando sus movimientos limitados a las direcciones contenidas 
en el plano tangente a la superficie por el punto de contacto. En tal caso el 
desplazamiento real no sería el natural que seguiría por acción del campo de 
fuerzas, sino el único compatible con el vínculo que lo liga a la superfície, 
pudiendo alcanzar el equilibrio no solo al anularse las fuerzas del campo sino 
también en caso de que dichas fuerzas fueran normales a la superficie, ya que 
z no habría componente en las únicas direcciones posibles de movimiento de las 
Figura 11.13 contenidas en el plano tangente. 
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Por el PROCEDIMIENTO DE LAS FUERZAS estas LIGADURAS encuentran su equivalente mecánico en unas 
REACCIONES que se originan para anular el movimiento en la dirección determinada (la normal a la superficie en 
el ejemplo), no para crearlo, y su valor Ñ es por tanto el necesario exclusivamente para anular la componente de 
las fuerzas activas en esa dirección. Las reacciones se consideran, por tanto, como FUERZAS PASIVAS "creadas": 
para sustituir mecánicamente las coacciones materiales, que impiden el movimiento pero no son causa por sí 
mismas, (de ahí su nombre en contraposición con las fuerzas activas de las que depende su existencia) y que se 
cuantifican de modo que su dirección viene fijada por la dirección del movimiento que impiden, su sentido es el 
contrario al que impiden y su módulo el necesario para lograr su objetivo "hasta alcanzar los límites de la 
deformabilidad". 


Según esta definición, las fuerzas pasivas pueden individualizarse sin problema para cada uno de los 
grados de libertad que limiten, siendo las más simples las que limitan un solo, posibilitando el movimiento en los 
otros dos (superficies), luego están las que limitan dos grados de libertad posibilitando únicamente el 
desplazamiento por una línea y, por último las que fijan la partícula a un punto limitando los tres grados de 
libertad. 


¿"Superficie (2 grados de libertad) Línea (1 grado de libertad) ' o (0 grados de libertad) 


Figura 11.14 


Ahora bien, podemos encontrar algunas diferencias según limiten unilateralmente el movimiento en un 
sentido pero no en el otro: LIGADURAS INCOMPLETAS (por ejemplo los apoyos), o lo impidan en cualquier sentido: 
LIGADURAS COMPLETAS. Matemáticamente se reflejará como una desigualdad en la función que representa la 
ligadura para las incompletas o en una igualdad para las completas: 


M4 


ey z)>0! 


I 
pxyz)<0 | os 


Ed h 7 Ligadura completa p(x yz)-0 
Ligadura incompleta p(x y z) < 0 . dl 


Figura 11.16 Figura 11.15 


En conclusión, las reacciones que intervienen en el planteamiento clásico de la Mecánica a partir de las 
Fuerzas se podrán expresar matemáticamente como vectores normales a las superficies y, por tanto, en la forma 
Ñ = hh grad y, si el punto está ligado a una superficie, o como dos vectores ÑN, =2, grad q, y N, = %, grád q, si 
el punto está ligado a una curva (intersección de dos superficies) y las ecuaciones desarrolladas para unos ejes 
cartesianos, teniendo en cuenta que se introduce una incógnita nueva, A¡por cada reacción N, a considerar, quedan: 
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Punto ligado a una Superficie: 


9% - 
F(xy2) +45 =0 


x=x 
a A F Z)+A2L=0 S (11.29) 
F(xyz)+2gradp=0 => y is dy 4 ecuaciones con 4 incógnitas > 3 e 
F,.(xyz)+1 22 =0 MG 
dz A=4 
p(xyz)=0 
Punto ligado a una Curva: 
op: 04. 
F.(xyz)+A e + Mor S 
99, 0 _ x=X 
) Ñ EJCIDA A hd (11.30) 
= = mx : .: A : e 
F(xyz) +2 grado, + A, grado, =0 > F(xy2) +2, ca +2, de a 5 ecuaciones con 5 incógnitas — Z=Z, 
p, (xyz)=0 Ay» A, 


p, (xyz)=0 


Sin embargo, podremos prescindir de las reacciones si planteamos el problema por procedimientos 
analíticos de Trabajo, ya que no producen trabajo para los desplazamientos compatibles con los vínculos por ser 
perpendiculares a ellos. El equivalente analítico de las ligaduras por este procedimiento son, en consecuencia, la 
propia expresión de la ligadura que limitará los desplazamientos a los únicos compatibles y nos permitirá 
relacionar las variables independientes de referencia para adecuar su número al número resultante de grados de 
libertad. 


Las ecuaciones equivalentes de la Estática para uno y otro método serán: 


F:9F=0 E, +hgrado(xyz) = 0 
Procedimiento de los Trabajos > Procedimiento de las Fuerzas 
p/(xyz)=0 0 (xyz)=0 


y como veremos a continuación, el Teorema de Lejeune-Dirichlet es válido también, si las fuerzas derivan de un 
potencial, para el punto ligado a una superficie o una curva, sin tener en cuenta la fuerzas de reacción, ya que la 
posición del punto y, por tanto, la variación del potencial, dependen exclusivamente de las fuerzas "activas" y de 
las limitaciones impuestas por el vínculo cuya ligadura consideraremos perfectamente lisa para no tener en cuenta 
tampoco las fuerzas debidas al rozamiento. 


PUNTO LIGADO A UNA SUPERFICIE D(x, y,z)=0 Si la resultante F' de las fuerzas aplicadas deriva de un potencial, 
como el gradiente de éste es normal a la superficie equipotencial correspondiente y, en los puntos de equilibrio, es 
también normal a la superficie D=0, solo habrá equilibrio en los puntos en que la superficie equipotencial que pase 
por ellos sea tangente a la superficie sobre la que el punto está obligado a permanecer. Estos puntos de tangencia 
corresponden, en general, a valores extremos de V, entre los que puede tomar esta función en los puntos de la 
superficie. Si V es máximo, el equilibrio será inestable, siendo estable para V mínimo, como definimos en el 
epígrafe anterior. 


PUNTO LIGADO A UNA CURVA O, (x,y,z) = 0, D,(x,y,z)= 0 También aquí es aplicable el teorema de Lejeune- 
Dirichlet en los mismos términos que para los casos anteriores del punto libre o ligado a una superficie. 
Análogamente a lo que ocurre para éste último, el equilibrio será únicamente posible si la tangente a la curva está 
contenida en un plano tangencial en P a la superfície equipotencial que pase por él. Tales puntos corresponden a 
máximos y mínimos del potencial V, entre los valores del mismo asequibles en el movimiento a lo largo de la curva. 
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